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ABSTRACT 


This report describes a Laser Doppler Velocimeter capable 
of measuring all three components of velocity simultaneously 
in low-speed flows. All the mean velocities, Reynolds stresses, 
and higher-order products can then be evaluated. The approach 
followed is to split one of the two colors used in a 2-D system, 
thus creating a third set of beams which is then focused in 
the flow from an off-axis direction. The third velocity component 
is computed from the known geometry of the system. In this 
report, the laser optical hardware and the data acquisition 
electronics are described in detail. In addition, full operating 
procedures and listings of the software (written in BASIC and 
ASSEMBLY languages) are also included. Some typical measurements 
obtained with this system in a vortex/mixing layer interaction 
are presented and compared directly to those obtained with a 
cross-wire system. A brief discription of the present system 
together with a review of existing 3-D Laser Doppler Velocimeters 
is given in Ref. 1. 
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NOMENCLATURE 


a(i) 

d f 

f 

1 

NS 

R 

S<a(i) > 

Sa 

u,v,w 

u,v,w 

u o 

x,y,z 

Subscripts: 

bragg 

mix 

dopp 

res 

1 

2 

Superscripts 


ith sample value 
fringe spacing 
frequency 
sample number 
number of samples 

range setting on Macrodyne processor 

sum of a(i) from i = 1 to i = Ns 

average of quantity a, defined by equations (3) 

mean velocity components in x,y,z directions, 
respectively 

instantaneous velocity components in x,y,z directions, 
respectively 

U X - U 2 , velocity difference between the two streams 
in the mixing layer 

Streamwise, normal, and spanwise coordinate directions, 
respectively 

Bragg 

mixing 

doppler 

resultant 

quantity measured in high-speed side of mixing 
layer 

quantity measured in low-speed side of mixing layer 

• 

fluctuating quantity, e.g., u = u + u' 
time average 
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INTRODUCTION 


Our ability to understand and model turbulent flows still 
relies heavily on the availability of accurate measurements 
of mean and fluctuating quantities within the flow. Until recently, 
the hot wire was the only reliable tool available for the measurement 
of fluctuating velocities. In fact, almost all of our present 
knowledge about turbulent flows is based on measurements made 
with hot wires. In relatively simple flows (moderately two- 
dimensional with small cross-flows) , reliable and accurate hot 
wire measurements are now possible with fully automated data 
acquisition and reduction systems which minimize errors due 
to drifts in calibrations. An example of such a system is given 
in Ref. 2. 

However, as we turn our attention towards more complex 
turbulent flows, a need for more sophisticated measurement techniques 
has become apparent. These complex flows include those with 
compressibility effects, strong three-dimensionality (with steep 
mean gradients) , flow reversals, and time-dependent behavior. 
Since about the mid-sixties, the most popular alternative tool 
for measuring mean and fluctuating velocities in turbulent flows 
has been the Laser Doppler Velocimeter (LDV) . 

The most popular LDV arrangement used for wind tunnel measure- 
ments is the dual beam or fringe method. In this method, one 
of the laser lines is split into two lines of equal intensity 
which are then focused through a lens so that they cross over 
at the focal point. The flow is seeded with small particles 
(typically less than 3-4 pm in diameter) which follow the fluid 
motion. As these particles pass through interferometric fringes 
created by the crossed laser beams, light is scattered off them 
which is received by a photodetector. The frequency of this 
scattered light, along with a knowledge of the fringe pattern 
formed by the laser beams, provide the means to calculate the 
velocity of the particle. The fringe method, especially in 
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the forward scatter, off-axis mode, generally offers the best 
signal-to-noise ratios and spatial resolution. 

Although LDV systems are somewhat complex and tiresome 
to set up, they have certain advantages over hot wires for turbulence 
measurements . The fact that Laser Doppler Velocimetry is non- 
intrusive is especially beneficial in the measurement of unstable 
flow phenomena which are very sensitive to the presence of measure- 
ment probes. In certain situations, LDV systems can also provide 
greater spatial resolution and better directional descrimination 
than hot wires. This makes it possible to use LDV systems for 
the measurement of separated flows. Since an LDV measures the 
velocity directly, independent of the thermodynamic properties 
of the flow, it is particularly attractive for velocity measurements 
in compressible flows. Furthermore, the calibration converting 
the frequency to velocity is linear and easy to implement in 
software. This feature also allows for uniform sensitivity 
in measuring both moderate and high turbulence intensities. 

Two-color LDV systems capable of measuring two components 
of velocity simultaneously are now being widely used. However, 
the main interest in the present investigation was to study 
three-dimensional interactions where it is desirable to obtain 
measurements of all three velocities. Hence, the first objective 
was to develop a laser velocimeter system capable of measuring 
all three components of velocity simultaneously so that all 
six components of Reynolds stress may be computed. Another 
objective was to compare these LDV measurements directly with 
those obtained with hot wires in flow fields where both techniqes 
are expected to perform satisfactorily. The purpose of this 
is to evaluate the performance of the system quantitatively 
and objectively. 

The approach followed in the present investigation was 
to convert an existing two-component LDV system into one capable 
of measuring all three components of velocity simultaneously. 
The system utilizes two wavelengths (488.0 and 514.5 nm) from 
a 4-watt Argon-Ion laser. The main four-beam matrix measures 
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u and v directly. The green line in the four-beam matrix is 
split (in half) using a dichroic filter and directed over the 
top of the traverse mechanism with mirrors, giving the third 
beam pair for the measurement of the w-component. This pair 
of beams (with rotated polarization) measures w sin 45°+ v cos 45°. 
Scattered light is collected in the off-axis forward scatter 
mode using two collection lenses. 

Signal processing is accomplished with single-particle 
burst counters, and the validated data are multiplexed through 
a "home-built" interface to an HP 9845B desk-top computer. 
Some selected first and second order products are reduced on-line, 
and the raw data is dumped onto floppy disk. An off-line program 
reduces the data, giving up to third order quantities and also 
plots histograms of the raw data for each channel. The software 
includes the capability to filter out noise by examining the 
histograms . 

The optical and signal processing hardware is described 
in Section 2. The data acquisition and reduction software is 
described in Section 3 and detailed operating procedures are 
given in Section 4. Some problems, inherent to 3-component 
LDV systems, are presented in Section 5. Sample results from 
an experiment measuring mean and turbulence quantities in a 
vortex/mixing layer interaction are compared directly to results 
obtained using crossed hot-wire anemometry in Section 6, and 
concluding remarks are presented in the final section. Complete 
software listings written in BASIC and ASSEMBLY languages to 
run on the HP 9845B desk-top computer are included in the appendix. 

2. OPTICAL SYSTEM AND SIGNAL PROCESSING HARDWARE 

The hardware for the 3 -component LDV system can be divided 
into three categories: the optical system, the signal processing 
instrumentation, and the computer. The LDV optics consist of 
the optics table, where the laser beam is split into green and 
blue beam pairs, the transmitting optics, where the beams are 
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directed into the flow field, and the receiving optics, where 
the scattered light is picked up by photodetectors. The signal 
processing instrumentation consists of amplifiers, filters, 
burst counters and a computer interface. The sampling procedures 
are all computer controlled. 

2.1 Optics Table 

The optics table consists of the laser and all the optical 
elements needed to provide the necessary four-beam matrix. 
Fig. 1 shows a schematic of the arrangement with component numbers 
as referred to in this section. A 4-watt Argon-Ion laser (Lexel 
Model 95) is used to produce the main beam. The beam then passes 
through a collimator (1) , which ensures that the beam waist 
occurs at the focal point of the transmitting lens . 

The collimated beam is passed through a color separator 
box, which consists of a polarization rotator (2) , an attenuator 
(3) , and a pair of high dispersion Brewster angle prisms (4) 
which are used to separate the multi-line beam into two colors, 
blue (488 nm) and green (514.5 nm) . These two beams are then 
reflected across the box by mirror (5) , and out of the box by 
mirrors (6, 7) . 

Following the color separation, the green beam’s polarity 
is rotated to horizontal (8) , and the beam is split into two 
beams in the vertical plane (9) . Most beam splitters prefer 
this type of perpendicular polarization for maximum efficiency. 
Using the beam displacer (10) , the blue beam is then moved to 
the center of the optics, and its polarity is rotated to vertical 
(11) . The blue beams are split in the horizontal plane (12) . 
At this point, the four beams are each displaced 25 mm from 
the optical axis. 

Two Bragg cells (13, 14) are used to shift the frequency 
of one beam from each pair. The unshifted beam passes through 
an optical rod so that the path lengths are matched. The frequency 
is shifted by a fixed amount of 4 0 Mhz . This shift creates 
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a moving system of fringes at the beam intersection point, allowing 
for directional discrimination of the velocity. Frequency shifting 
also helps to reduce the percentage of frequency change in highly 
turbulent flow, to reduce fringe bias, and to optimize frequency, 
thus enabling easy removal of the pedestal by high pass filtering. 

The four beams are then passed through a beam steering 
module (15) . The module consists of a set of wedge prisms that 
can be independently rotated about the beam axis to steer the 
shifted beam in any direction. This allows a more precise alignment 
of the beams. The beams are finally passed through a beam displacer 

(16) to reduce the beam spacing to 13 mm and a rotating prism 

(17) before leaving the optics table;. The rotating prism enables 
the four-beam matrix to be rotated independently so that the 
beams may be aligned relative to the tunnel axes. 

2 . 2 Transmitting Optics 

The transmitting optics (Fig. 2) are mounted on a traversing 
mechanism with three degrees of freedom. The traverses are 
driven by individual stepper motors. The four beams from the 
optics table are directed by a set of five mirrors through a 
dichroic filter before being focused by a 380 mm (15 inch) focal 
length lens. This main four-beam matrix measures u and v directly. 
The dichroic filter, set at an angle of about 30 degrees to 
the incoming beams, splits the green beams in half, which provides 
the third beam pair for the third velocity component, w. This 
third beam pair is directed over the top by mirrors, passed 
through a polarization rotator (giving it a different polarity 
than the main-axis green beam pair) and then focused at the 
focal point of the main beam set. Since this third beam pair 
intersects the main measuring volume at a 45° angle to the main 
axis, it measures v and w with equal sensitivity, with the measured 
component being w v = (v cos 45° + w sin 45°) . Since v is measured 
directly, w can be evaluated using the equation: 
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w 


(1) 


w v - v cos 45° 
sin 45° 

Typical probe volume dimensions for each beam pair in the 
present configuration are 10 mm in length and 0.2 mm in diameter 
(Fig. 3) . However, the actual "viewed" dimensions are reduced 
considerably, as discussed below, in Section 2.3. 

2 . 3 Receiving Optics 

The detector system (Fig. 4) is in the off-axis forward 
scatter mode. The receiving optics are mounted on a traversing 
gear, also run by stepper motors, which moves synchronously 
with the transmitting optics. Scattered light is collected 
by two 380 mm (15 inch) focal length lenses. The collimated 
light is passed through filters to separate the colors and the 
off-axis line is additionally passed through a polarization 
filter to avoid collecting light scattered by the main axis 
green pair. The collected light is then focused by 250 mm (10 
inch) focal length lenses onto pin-hole apertures mounted in 
front of the three photomult ip ler tubes. The collection angle 
and diameter can be adjusted to select the effective (viewed) 
probe length (Fig. 5) . In the present set-up, a collection 
angle of 3 0 degrees and an aperture diameter of 0.5 mm were 
used to give an effective length of about 1.5 mm. 

2.4 Signal Processing Hardware 

The signals from the photomultiplier tubes are amplified 
and relayed to the signal processors via high-pass filters and 
mixers (Fig. 6) . The amplifiers used are EIN model 403LA with 
a fixed gain of 37 dB and the filters are Allen Avionics F2440 
with a fixed high-pass cutoff of 10 MHz. The filtered signals 
are mixed electronically with sine waves from three Tektronix 
SG503 Levelled Sine Wave Generators. The mixer is commercially 
available from Hewlett-Packard, model 1053 4A. The mixing procedure 
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is necessary for low-speed flows, where the actual Doppler fre- 
quencies are small compared to the Bragg frequency of 40 MHz. 
So in order to reduce the effective measured frequency and hence 
improve the counter resolution, the incoming signals are mixed 
with sine waves of known frequency. 

The mixed signal is fed into single particle burst counters 
(Macrodyne model 2096-2 and 3003) via high-pass/low-pass filters 
and an amplifier (xlO) . The counters measure the zero crossings 
of the Doppler signal, which is related to the Doppler frequency 
by the range set on the Macrodyne. The range is set manually 
based on expected flow velocities, since it limits the frequencies 
that the processor can see for the given 10 bits of resolution. 

The processors use two checks to validate a Doppler signal. 
The first check is the usual 5/8 comparison, where the processor 
checks the frequency for 5 zero crossings against that for 8 
crossings. The second check is the multi-sequence check. Positive 
and negative thresholds are set on the signal, and a validated 
output is permitted only if, for all eight fringe crossings, 
the ■ signal passes through a positive threshold, a zero level, 
and a negative threshold in the proper sequence. The digital 
data (consisting of a 10 bit data word with 3 bits giving the 
range) and a sync pulse (produced every time the front end of 
a valid burst is detected) are passed to the computer interface. 

2 . 5 Computer Interface 

A NASA LDV-A/D computer interface (Cl) is used to transfer 
data from the LDV signal processor to the computer (Fig. 7) . 
The Cl can interface either digital or mixed analog and digital 
data to an HP 9 8 45B desk-top computer. The Cl consists of an 
eight-channel multiplexer, a four-channel A/D converter, and 
an event synchronizer with time interval counter. 

For use with the LDV, the inputs to the Cl are all digital. 
Six of the eight words come from the processor, and two are 
time and status words from the synchronizer. The inputs are 
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multiplexed to a single digital data channel output. 

The Cl can accept data from the three processors in either 
random mode or sync mode. In random mode, the Cl will accept 
data inputs when an event occurs on any of the three channels. 
A dead time (between 5 and 50 ps) is set in this mode, which 
controls the minimum time between samples to ensure that a given 
particle is sampled only once. In sync mode, the Cl will accept 
data inputs when all three processors sample simultaneous events. 
In this mode, a coincidence time (between 5 and 50 ps) has to 
be set, which determines the time window within which all three 
events must occur. A detailed description of the Cl can be 
found in Ref. 5. 

2.6 HP 9845B Desk-top Computer 

Multiplexed data are passed to the HP computer from the 
NASA LDV-A/D using the HP 98032A high-speed 16-bit parallel 
interface. Jumpers labeled "9,B,D" are connected inside the 
98032A for proper operation with the LDV-A/D. A select code 
of 10 (screw setting on the 98032A) is set for use with the 
software described in the appendix. A data buffer of 24 kbytes 
is provided in the memory of the HP 9845B desk-top computer 
for storage of up to 2000 samples obtained from one measurement 
location. Since three of the six data words passed for each 
sample are merely monitored and discarded as described above, 
12 kbytes of raw data remain to be stored for each point. An 
HP 9895A floppy-disk drive is used for archival storage of raw 
data. The buffered data are written in real time to a sequential- 
access floppy-disk file. Enough header information is written 
to each file to identify the run, and to reproduce calibration 
tables. 

3. DATA ACQUISITION AND REDUCTION SOFTWARE 

Data aquisition and reduction on the HP 9845B is done via 


8 



two programs. The data acquisition program controls the processor 
sampling. The program accepts heading and initialization parameters 
provided by the operator, performs a fast I/O handshake to acquire 
the raw data from the LDV-A/D computer interface, and writes 
the raw data and initialization parameters onto floppy-disk. 
In addition, the data acquisition program has a limited capability 
for on-line data reduction so that key results may be monitored 
during a run. The off-line data reduction program converts 
the raw data into instantaneous velocities and computes all 
the statistical quantities. The off-line program also has options 
to plot histograms of the raw data, to plot profiles of the 
reduced data, and to filter out noise. 

3.1 Data Acquisition Program 

The data acquisition program structure can be divided into 
four areas: (1) initialization of variables, (2) data acquisition, 
(3) limited data reduction, and (4) raw data storage (see the 
block diagram in Fig. 8) . During the initialization stage of 
the program, relevant test parameters are input to be saved 
along with the data. These parameters include the spatial coord- 
inates, Bragg and mixing frequencies, and angle of the third 
beam pair. Fringe spacings and probe dimensions are also calcu- 
lated at this time. 

Data is acquired from the LDV computer interface by means 
of a fast I/O handshake. The operator requests that a certain 
number of data samples be taken and stored in the data buffer 
as described in Section 2.6. The 16-bit data words are put 
into an integer array, three words wide, containing values for 
the three velocity components. The number acquired from the 
counter occupies only bits 0-9 of the data word. Bits 10-13 
specify the range, and bits 14 and 15 are unused. The data 
array, along with the initialization parameters given by the 
operator, are stored in a floppy-disk file. As an option, some 
of the data can be reduced on-line, before the raw data is written 
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to floppy-disk. This on-line reduction feature uses the same 
algorithms as the off-line data reduction program, but is faster 
since it computes fewer quantitites. Data acquisition times 
depend on many factors other than the software, but on-line 
data reduction takes about 30 seconds per point on the HP 9845B, 
and each file write to the floppy-disk requires about 20 seconds, 
assuming 2000 samples are taken per point. 

3.2 Data Reduction Program 

A separate off-line program is used to reduce the complete 
data. The data reduction program has four major sections: 
(1) data read from floppy-disk, (2) conversion of raw data to 
instantaneous velocities, (3) calculation of sums of instantaneous 
velocity components, and (4) calculation of moments from the 
sums (see the block diagram in Fig. 9) . Since the HP 9845B 
is a relatively slow micro-computer, the bulk of the calculation 
routines are written in ASSEMBLY language to reduce computation 
time. 

The program reads the initialization parameters and raw 
data which were written on floppy-disk by the data acquisition 
program. The raw data is used together with the mixing frequencies, 
Bragg shift frequency, third beam crossing angle, and fringe 
spacings to calculate the three velocity components for each 
sample. . The velocities for each sample are calculated in turn, 
and a running sum is maintained. After the samples have been 
summed, the average values of the various moments are computed, 
and the mean velocities and turbulence quantities are calculated 
from them. Reading from the floppy-disk takes about 20 seconds, 
while the data reduction including histogram plots requires 
about 90 seconds per point on the HP 9845B. 

The procedure for data reduction is fairly straightforward. 
The raw data is converted to frequencies using the range set 
on the processor. The measured frequency is given by the relation 
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= 3.2 * 10 4 / 2 R * D 


( 2 ) 


^res 

where D = raw data. 

The velocities are then calculated from the resultant frequen 
cies, by the relation: 


velocity = d f * ((fbragg- f mix) “ f res) ( 3 ) 

Once the velocities have been evaluated, average values of the 
various moments are computed as defined below: 


Su = S<u(i)> / Ns 
Sv = S<v(i) > / Ns 
Sw = S<w(i)> / Ns 
Suu = S<u(i)u(i)> / Ns 
Suv = S<u(i) v(i) > / Ns 
Suw = S<u(i)w(i)> / Ns 
Sw = S<v(i)v(i)> / Ns 
Svw = S<v(i)w(i)> / Ns 
Sww = S<w(i)w(i)> / Ns 
Suuv = S<u(i)u(i)v(i) > / Ns 
Suuw = S<u(i)u(i) w(i) > / Ns 
Suw = S<u(i) v(i) v(i) > / Ns 
Suww = S<u(i)w(i) w(i) > / Ns 
Suvw = S<u(i) v(i) w(i) > / Ns 


( 4 ) 


Using these definitions, the signal statistics are then calculated 
assuming nearly infinite sample size: 


u 

= 

su 


V 

= 

Sv 


w 

= ’ 

SW 


u' 2 

= 

Suu 

- SuSU 

V 2 

— 

Sw 

- SvSv 
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Sww - SwSw 


w 


,2 


u'v' 


u'w' 


v'w' 

u' 2 v' 

u ,2 w' 

u'v' 2 

u'w* 2 


U« v'w' 


Suv - SuSv 
Suw - SuSw 
Svw - SvSw 

Suuv - 2SuSuv - SvSUU + 
Suuw - 2SUSUW - SwSuu + 
Suw - 2SvSuv - SuSw + 
Suww - 2SWSUW - SuSww + 
Suvw - SuSvw - SvSuw - 


2SvSuSu 

2SwSuSu 

2SUSVSV 

2SuSwSw 

SwSuv + 2SuSvSw 


The implementation of the data reduction software is somewhat 
more complex. In order to achieve reasonable running times, 
the ASSEMBLY code is optimized for speed rather than clarity 
of operation. Operations which do not change between samples 


are performed only once. Thus, 

velocity = d f * ((fbragg " fmix) " f res) * 106 ( 6 ) 

becomes 

f int = f bragg “ f mix > d int = df * 10 6 (7) 

hence, velocity = dj. nt * (fi n t - f res) ( 8 ) 

Additional calculations are required to obtain the w component 
velocity. Since the third set of beams measures w v = w sin 45° 
+ v cos 45° , to find w, we must also perform the calculation: 

w = (w v - v cos 45°) / sin 45°. (9) 


The raw data output of the digitizer is in the form of 
a 10 bit data word, giving a range of possible values (counts) 
from 0 to 1023. Since 2000 samples are normally taken for each 
point, a given value may be encountered many times. Accordingly, 
each time a new value of the raw data is encountered, the cor- 
responding velocity is calculated and stored in a look-up table. 
The next time that value is encountered by the program, the 
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proper velocity can be easily looked up, eliminating the need 
for another time-consuming real variable calculation. The exact 
running procedure of the data reduction program therefore, is 
as follows: 

1. Reads header and all raw data for a particular point from 
floppy disk. 

2 . Strips bits 10-13 from three raw data words corresponding 
to the three different velocity components, and uses them to 
calculate the three ranges set on the A/D. 

3. Calculates various intermediate values which remain the 
same throughout the point. 

4. Reads a data word, strips off bits 0-9, and uses this raw 
datum as an index to look up its corresponding velocity. 

5. If it is a new value, the program calculates the velocity 
using the raw datum, and stores it in the appropriate place 
in the look-up table. 

6. For the w velocity component, it finds the actual velocity 
from w = (w v - v cos 45° ) / sin 45° . 

7. Performs steps 4, 5, and 6 three times — once for each velocity 
component of the sample. 

8. Updates the running sums of the velocity components and 
products of the velocity components. 

9. Performs steps 4 through 8 for all samples. 

10. Uses the sums to obtain the average velocities, Reynolds 
stresses, and third order products, and then prints these quantities. 
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11. Plots histograms of the raw data for all three channels. 

12. Performs steps 1 through 11 for each profile point. 

13. Plots profiles of the reduced data. 

14. Tabulates normalized data profiles. 

15. Writes a summary file containing the reduced data to disk. 


The data reduction program also has a routine to filter 
noise and spurious data. Each of the 2000 data samples has 
3 counts corresponding to the 3 velocities associated with it. 
The filtering routine causes the data reduction program to ignore 
samples associated with counts which are excessively far from 
the mean. 

The filtering routine first sorts the data into three frequency 
tables, one for each channel. In a table, each count, i, has 
associated with it a number, Sj_, which is the number of samples 
with that particular count. 

The filtering routine finds the average count for each 
channel by going through the frequency tables and using the 
formula 

1024 

T . 

Cj — i=i 1 — x (10) 

Ns 


where j = 1,2,3 = channel number 


Cj = average count 
Ns = number of samples. 

Next, the filtering routine finds the standard deviation by 


performing the summation 

1024 





(ID 
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where 0" j = standard deviation. 

An input variable called S^ ev is read by the filtering 
routine. All counts further than S^ ev standard deviations away 
from the mean will be filtered out. This is done by multiplying 
s dev b Y CJj for each channel , and going through the frequency 
tables one more time. If the magnitude of (i - Cj) > S^ev * CT j / 
then S-^j is set equal to zero. 

The routine which calculates the statistical quantities 
takes each sample one at a time. The three counts associated 
with each sample are found and looked up in the frequency table. 
If for any count the table entry (Sj_j) is zero, then the sample 
is discarded. Thus any sample which is excessively far from 
the mean in any one of its three counts is not used. 

4. OPERATING PROCEDURE 

4.1 Alignment Procedure 

To achieve the best possible beam crossing and the most 
effective measuring volume, each module in the optics system 
must be carefully aligned. Detailed alignment instructions 
for individual optical components are given in Ref. 3. The 
overall alignment procedure is described in this report. Component 
numbers in this section refer to those shown in Fig. 1. 

The first step is to check that the laser output beam is 
parallel to the optics table at the specified height, using 
the system alignment blocks. The collimator (1) should be positioned 
so that the laser beam goes through the center of the lens, 
and focused so that the beam waist occurs at the cross-over 
point. The aligned beam then passes through the polarization 
rotator (2) and into the color separator box. The attenuator 
(3) ensures that the beam has horizontal polarity at this point. 

Each component of the color separator must be aligned separ- 
ately. The dispersion prism (4) should be aligned so that the 
path length is the same in both prisms and is parallel to the 
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base of each prism. Several beams emerge from the prism, with 
the two brightest ones being green (514.5 nm) and blue (488 
nm) . The two beams are reflected out of the box using mirrors. 
The beams should be centered on mirrors (6) and (7) . As the 
beams come out of the box, alignment blocks are used to check 
the beam positions. If both beams are off-axis in the same 
direction, mirror (5) is used to align them. If only one beam 
is off, the appropriate mirror, (6) or (7) , is used for the 
adjustment. Beam splitter efficiency is maximized when the 
beam polarization is perpendicular to the plane of the split 
beams. This is achieved through components (8) and (11). The 
beam splitters (9) and (12) cannot split the beam intensity 
exactly in half, so one beam of each pair is always slightly 
brighter. Since the Bragg cells (13) and (14) normally attenuate 
the shifted beam, the brightest beam is frequency shifted so 
that the output beam pairs have nearly equal intensities. Detailed 
instructions for aligning the Bragg cells are given in Ref. 4. 
After passing through the Bragg cells, directional wedges are 
used to project the beams onto a distant surface (~ 3 m) . Any 
beam misalignments are more easily seen this way, and with the 
use of a marked mask, the beams can be adjusted to the correct 
orientation. Mirrors (6, 7) are used to adjust the unshifted 
beams, and the beam steering modules (15) are used for rotating 
the shifted beams along two circular arcs. 

The next check is to ensure that the beams are parallel 
to each leg of the traverse mechanism on the transmitting optics 
table, so that beam alignment is maintained while traversing. 
The dichroic filter (Fig. 2) is adjusted so that the green beam 
pairs are split equally. Once the beams pass through the trans- 
mitting optics, the four beams must be arranged so that they 
all cross at the same position. A microscope objective is used 
to view the beam crossing. If the four beams are not symmetric 
about the optical axis, mirrors (6) or (7) can be adjusted to 
correct this. The beam steering modules (15) are used to ensure 
that the beams cross at the same point along the axis. The 
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third pair of beams are now aligned (by eye) so that they also 
cross at the same point as the main line beams. 

To align the receiving optics, a piece of translucent tape 
is placed at the beam intersection point to scatter the laser 
light. By tracking the scattered light, the receiving optics 
are aligned to give a sharp image at the pin-hole aperature 
in front of the photomultiplier tubes. 

4.2 Signal Processing 

The sensitivity of the photomultiplier (PM) tubes used 
in the receiving optics can be varied by varying the applied 
voltage. Typically, a voltage of about 1000 volts is applied. 
This voltage can be increased to make the PM tubes more sensitive, 
as long as the threshold levels on the processors are increased 
accordingly, since the amount of noise picked up is also in- 
creased. 

The measured signal is mixed electronically with the signal 
from a sine wave generator. The frequency of the sine wave 
(the mixing frequency) is chosen based on the expected flow 
velocities. The mixing frequency is chosen such that the difference 
between it and the Bragg shift is about twice the maximum expected 
Doppler frequency. This allows enough margin for fluctuations 
about the expected Doppler frequency and still have a remaining 
nonzero resultant frequency (f rea = fbragg - f mix " f dopp) • 
If f mix i s to ° high, a biasing results, similar to the fringe 
biasing caused by stationary fringes. The number of fringes 
crossed by a particle per second (as seen by the processor) 
is proportional to At(fb ra gg - fmix)' where At is the time taken 
by a particle to cross one fringe. Noting that At is only determined 
by the fringe spacing and the flow speed, if f m i x is increased f 
the number of fringes crossed by a particle is effectively reduced. 
This means that signals from particles which cross the fringes 
at an angle may not have enough fringe crossings (8) to be validated, 
and hence a bias towards particles moving perpendicularly to 
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the fringes (higher velocity) results. This gives a higher 
mean velocity but a lower fluctuation level. 

The amplitude of the sine wave must be chosen so that an 
adequate signal-to-noise ratio is maintained. Typically, a 
peak-to-peak amplitude of 1 volt is required. 

The mixed signal is then fed into single particle burst 
(Macrodyne) processors. The high pass and low pass filter frequen- 
cies are set so that the processor frequency is centered between 
the two. Usually the high pass filter is set at 0.5 Mhz, which 
is the lowest non-zero setting, and the low pass filter is set 
anywhere from 2 to 16 Mhz, depending on the magnitude of the 
velocity component being measured. The filter bandwidth should 
be broad enough so that no parts of the fluctuating signal are 
attenuated. 

The gain on the processor is normally set to 10. The output 
of the processors is displayed on an oscilloscope, and the Doppler 
signals should read about 1 volt peak-to-peak. The PM tube 
voltage can be adjusted so that the signal is at the desired 
level. Signal levels of more than about 1 volt end up being 
clipped and will therefore not be validated by the processor. 

The comparator accuracy for the 5/8 signal validation test 
can be set between 0 and 10 count variation. The 0 setting 
is the most accurate, and 11 is off. The processor manufacturer 
recommends that this level be usually set to 9 . 

The range on the processor sets the bandwidth of frequencies 
that the processor can see, according to Table 1. For each 
range, a particular frequency corresponds to a count, from o 
to 1023. The processor frequency should be matched to a number 
in the central column in Table 1, and the corresponding range 
should be set on the Macrodyne processor. An additional check 
can be made on the range setting by monitoring the analog output 
from the Macrodyne on a DC voltmeter. The output ranges from 
0 to 10 volts, which corresponds to the counts from 0 to 1023. 
The correct range setting is that which gives about 5 volts 
on the voltmeter at the operating velocity. 
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The threshold is then set for the multi-sequence check. 
The threshold should be set so that the data rate is about half 
of the data rate at zero threshold. Another check is to block 
one beam in each color pair. The data rate on the corresponding 
channel should be zero when one beam is blocked. 

To obtain good data rates in air, the flow must be seeded 
with uniformly sized particles. Smoke, obtained from burning 
mineral oil or incense, provides particles of approximately 
2 pm. 

The computer interface should be set up as shown in Fig. 10. 
First choose the coincidence mode. If shear stresses are to 
be evaluated, coincidence on all 3 channels must be selected. 
The coincidence time (5 to 50 jis) sets the window width during 
which coincidence is defined. This should be set at 5 jis . 
If only the individual velocities are desired, then the random 
mode can be selected. The dead time (5 to 50 p.s) should then 
be set so that data from one particle is not recorded twice. 
A setting of about 25 jis is recommended for low-speed flows. 

The number of words that must be multiplexed can be calculated 
as follows: #words = #inputs + 2. This number should be rounded 
to the nearest even number. In the present case a setting of 
six is used. The event mode is set to LDV (digital data only) . 
The LDV- A/D switch enables both digital and analog data to be 
interfaced simultaneously. The counter clock frequency is set 
to equal the approximate data rate, and the computer select 
is set to HP. 

5 . INHERENT PROBLEMS 

Some design problems inherent to 3-component LDV systems 
are discussed in this section. One main problem with some of 
the earlier designs which called for splitting a color to create 
the third beam pair had to do with cross-talk. This is where 
signals from the two channels bearing the same color could not 
be adequately separated. In the original design of the present 
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set-up, the polarization of one of the green beam pairs was 
rotated relative to the other, so that the receiving optics 
could distinguish between the two signals. However, using a 
relatively large angle for the off-axis beams (45° ) and two 
separate collection lenses, cross-talk between the two green 
channels has been almost eliminated, thus making the polarization 
rotation dispensable. 

Some earlier designs of 3 -component LDV systems measured 
the u + w velocity component with the off-axis third beam pair. 
It is shown in Ref. 1 that measuring the v + w component instead, 
as done in the present system, reduces the uncertainty in the 
w component relative to the uncertainty in w from these earlier 
systems . 

Another problem has to do with signal coincidence. Details 
of the probe volumes for the present system are shown in Fig. 3. 
It is clearly illustrated how the cross-over region between 
the three sets of beams forms a very small fraction of the overall 
probe volume. Thus, with heavy seeding (necessary for three-channel 
work) , the electronics may validate data received from different 
particles which are not necessarily in the cross-over region 
but are within the coincidence time window set on the interface. 
This results in a lack of correlation between the measured veloc- 
ities, and causes the evaluated shear stresses to be inaccurate. 
Two schemes have been used in an attempt to minimize this problem. 
First, the coincidence time was made so short that measurements 
from different particles may be considered instantaneous. The 
minimum setting of 5 ;as available in the present hardware was 
used; this is equivalent to about half the flight time of a 
particle passing through the probe volume. The second procedure 
involved reducing the effective "viewed" probe length and thereby 
reducing the probability of this "apparent coincidence" (Fig. 5) . 
A collection angle 0 of 30° and an aperture diameter of 0.5 mm 
were used to give an effective length of about 1.5 mm. 
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6 


SAMPLE RESULTS 


As an initial check on the accuracy of the present system, 
some preliminary measurements have been made in a vortex/mixing 
layer interaction, previously investigated using the cross-wire 
technique (Ref. 6). Since the induced cross-flow angles in 
this interaction are only 5 P -10° , cross-wire measurements are 
expected to be accurate to within about 5%. A schematic of 
the experimental set-up is shown in Figs. 11a and b. LDV measure- 
ments of the secondary flow velocities at one streamwise station 
(x = 229 mm) are presented and compared to the cross-wire measure- 
ments in Figs. 12a and b. LDV and cross-wire measurements of 
the turbulence quantities at one spanwise position (z = 13 mm) 
are compared in Fig. 13. 

The secondary flow velocities are qualitatively similar, 

although the LDV measurements indicate a somewhat higher w. 

The normal intensity w' 2 also seems slightly high. The higher 

w measurements are more likely caused by a slight misalignment 

of the beams relative to the tunnel axis rather than by remnants 

of the apparent coincidence problem (discussed above in Section 5) , 

since the latter would not affect the w measurements, v' 2 seems 

to agree very well whereas u' 2 is a bit low, and since u was 

a bit high, this was probably a result of fringe biasing caused 

by too high a mixing frequency (as discussed above in Section 4.2) . 

(The mixing frequencies used for these measurements were 37.5 Mhz 

for the u channel, and 38 MHz for the v and w v channels.) However, 

the normal stress measurements agree to within 10%, and the 

shear stresses are consequently affected. u'v' is somewhat 

low (about 20%) whereas u'w' is low by almost a factor of two. 

The measurement of u'w' with the present system seems to be 

very sensitive. This is due to the fact the w v * u is of the 

o , 

same order as u * v cos 45 , so any small error in the measurement 
of these velocities or the off-axis angle can result in large 
errors in u'w'. v'w' in vortex affected flows is generally 
of the same order as u'w' and this seems to be the case with 
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the present measurements. With a cross-wire, v'w' has to be 
evaluated from measurements made in four different planes about 
the probe axis, and hence was not measured here. The comparisons 
clearly demonstrate the potential of the system in measuring 
detailed mean flow and turbulence quantities in three-dimensional 
flows. Work is in progress on optimizing the problems discussed 
above so that the measurement accuracies may be improved. 

6. CONCLUDING REMARKS 

A 3 -component LDV system, capable of measuring all three 
components of velocity simultaneously has been developed for 
use in low-speed three-dimensional flows. All the six components 
of Reynolds shear stress and higher order products of interest 
can hence be evaluated. The approach followed was to convert 
an existing 2-component system by splitting one of the colors 
to produce the third beam pair. The additional optical hardware 
required for this process is relatively minor. 

For the first time, three- component measurements made with 
an LDV system have been compared directly with those obtained 
with the cross-wire technique, in a three-dimensional flow field 
where both techniques are expected to perform satisfactorily. 
The preliminary measurements are encouraging and work is in 
progress on improving the system accuracy. 
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APPENDIX 


SOFTWARE FOR THE HP 9845B DESK-TOP COMPUTER 

Complete listings of two programs written in BASIC and 
ASSEMBLY languages are included in this appendix: "LDV" for 
data acquisition, some on-line data reduction, and storage of 
data on floppy disks; and "STAT" for complete off-line data 
reduction from files written to disk. 
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REM PROGRAM LD7 

! PROGRAM TO ACQUIRE DATA FROM THREE-COMPONENT LBV SYSTEM 
! The program asks for initialization data and calculates the 
! calibration constants from them. It reads 3 channels of raw 
! LV data from the LBV A/D Cl and writes them to a disk file 
! together with the calibration constants. The program can reduce 
! the data and display real-time histograms of the raw LV data 
! if desired. 

OPTION BASE 1 
ICOM 16600 
IDELETE ALL 
I ASSEMBLE Find_uel 
I ASSEMBLE Data trans 
I ASSEMBLE Drawjiist 

COM INTEGER Dat a< 3 , 2000 > , D 1 < 2000 , 6 > , Ns , Nn ! D1 is data buffer 
COM REAL Df 1 , Df 2 , Df 3 

COM REAL Theta, Nub, Numi xl , Numi x2, Numi x3 
COM REAL Su, So, Sw, Suu, Sou, Sww, Suo, Suw, Sow 
DIM Date* C80] , F i 1 e*[6] , Name*[4] , Ti t 1 *[80] 

REAL Xpos , Ypos , Zpos 
INTEGER A, B, Rangel , Range2 
INTEGER Run, Dn, Mss, N 
REAL Phi , D, F, Ph2 

REAL Laml , Lam2, Lam3, Db, Prwi dl , Pr 1 enl 
REAL Re, Ue 

REAL Prwi d2, Prw i d3 , Pr 1 en2 , Pr 1 en3 

REAL Nfr 1 , Nfr2, Nfr3 

REAL U,V,W 

Ns=2000 

DEG 

PRINT 

PRINT " ** << PROGRAM LDV : 3-COMPONENT VELOCITY DATA >> ** " 

PRINT 

PRINT "PROGRAM STRUCTURE" 

PRINT " 1. INITIALIZE VARIABLES AND CALCULATE PARAMETERS" 

PRINT' " 2. ACQUIRE DATA FROM A/D" 

PRINT " 3. WRITE TO FLOPPY DISC" 

PRINT 

! 

! ** CHECK HISTOGRAMS ** 

! 

Ans*= " N " 

INPUT "DO YOU WISH TO LOOK AT HISTOGRAMS ? <Y/N, DEFAULT N)",Ans* 
IF Ans*="Y" THEN GOSUB Hist 
! 

! ** INITIALIZE RUN ** 

! 

PRINT " ** INITIALIZATION ** " 

Run=l 

PRINT " ENTER RIJN PARAMETERS:" 

PRINT 

INPUT "Enter date and time:", Date* 

INPUT "Enter 1-Line Name For Profile : " , T i 1 1 * 

INPUT "No. of data samples per point <2000 maxm. > : " , Ns 

I 

! BEAM SPACING D IS FIXED 
! 

D= . 0 1 3 

j 

! FOCAL LENGTH F IS FIXED 
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610 ! 

620 F= . 38 1 
630 ! 

640 ! WAVELENGTHS OF 3 BEAMS ARE FIXED 

650 ! 

660 Laml=4. 83E-7 
670 Lam2=5. 145E-7 

680 Lam3=5. 145E-7 

690 ! 

700 ! REFERENCE VELOCITY SET TO ZERO 

710 ! 

720 Ue=0 
730 ! 

740 ! CALCULATE HALF-ANGLES FROM BEAM SPACINGS RND FOCAL LENGTH 

750 ! 

760 Phl=ATN<D/'2/F> 

770 Ph2=Phl*2 
780 ! 

790 ! CALCULATE FRINGE SPACINGS FROM WAVELENGTHS AND HALF-ANGLES 

800 ! 

810 Df l=Laml/'2/SIN<Phl) 

820 Df2=Lam2/2-'SIN<Phl> 

830 Df3=Lam3/2xSIN<Phl) 

840 ! 

850 ! BEAM DIAMETER Db IS FIXED 

860 ! 

870 Db= 1 . 20E-3 
880 ! 

890 ! CALCULATE PROBE VOLUME WIDTH AND LENGTH 

900 ! 

910 Prwidl=4*Laml*F.'PI/Db/'C0S<Phl> 

920 Prlenl=4*Laml*F-'PI/Db/'SIN<Phl> 

930 Prwid2=4*Lam2*F.'PI/'Db/'C0S<Phn 
940 Prl en2=4*Lam2*F.'PI/Db/'SIN<Phl> 

950 Prwi d3=Prui d2 
960 Pr 1 en3=Prl en2 
970 Nfrl=Prwidl/Df 1 
980 Nfr2=Prui d2xDf2 
990 Nfr3=Nfr2 
1000 ! 

1010 ! GET MORE RUN PARAMETERS 

1020 ! 

1030 INPUT "Enter Bragg shift frequency <MHz) : " , Nub 
1040 Ki 1 l_w=3 

1050 INPUT "Is this a two-channel or three-channel run (2/3, default 3) ?",Kill 
w 


1060 

IF Ki 1 l_w=3 THEN GOTO 1110 



1070 

INPUT "Enter mixing frequency (MHz, 

2 

nos . > 11 , Numi x 1 , Num i x2 

1080 

Numi x3=0 



1090 

Theta=0 



1100 

GOTO 1130 



1110 

INPUT "Enter mixing frequency (MHz, 

3 

nos. ) " , Numi xl , Numi x2, Numi x3 

1120 

INPUT "Enter 3rd beam angle (degrs) 

• II 

, Theta 

1 130 

INPUT "Enter tunnel reference voltage 

< volts)" , Vref 

1140 

PRINTER IS 0 



1 150 

j 



1160 

! PRINT HEADER 



1170 

! 



1180 

PRINT Ti 1 1 $ 



1190 

PRINT "TEST DATE AND TIME :",Date$ 



1200 

PRINT "BEAM SPACINGS (m> = ",D 
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1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1599 

1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 


PRINT "TOTAL ANGLE BETWEEN BEAMS CDegrs) =" ,Ph2 
PRINT "FRINGE SPACINGS Cm) = " , Bf 1 , Df 2 , Df 3 
PRINT "PROBE WIDTHS Cm) = " , Prw i d 1 , Prw i d2 , Prw i d3 
PRINT "PROBE LENGTHS Cm) = " , Pr 1 en 1 , Pr 1 en2 , Pr 1 en3 
PRINT "NO. OF FRINGES = " , Nfr 1 , Nfr2, Nfr3 
PRINT "BRAGG SHIFT FREQUENCY CMHz) =",Nub 
PRINT "MIXING FREQUENCY CMHz) =" , Numi xl , Numi x2, Nurni x3 
PRINT "THIRD BEAM SET ANGLE CDegrs) =", Theta 
PRINT "RUN NUMBER = " , Run 

PRINT "TUNNEL REFERENCE VOLTAGE Cvolts) =" ,Vref 

PRINTER IS 16 

Ans$="N" 

INPUT "DO YOU WISH TO MAKE ANY CHANGES? CY^N, DEFAULT N)",Ans$ 
IF Ans$="Y" THEN GOTO 540 
! 

! DATA ACQUISITION 
! 

PRINT 

PRINT " ** DATA ACQUISITION ** " 

REDIM DataC3,Ns> 

Ans$="N" 

INPUT "DO YOU WISH TO CHANGE FILE NAME ? CY/N, DEFAULT N)",Ans$ 

IF Ans$="N" THEN GOTO 1490 

INPUT "Enter 4-digit f i 1 enamel ", Name* 

PRINTER IS 0 

PRINT "4-digit filename for profile :",Name$ 

INPUT "Enter Disk Number : " , Dn 
PRINT "Disk Number =",Dn 
PRINTER IS 16 

INPUT "Enter X, Y, AND Z 1 ocat i ons: " , Xpos, Ypos, Zpos 
! 

! TAKE DATA 

j 

PRINTER IS 0 

PRINT 

PRINT 

PRINT "POINT NUMBER IN PROFILE : ",Run 
PRINT " X , Y , Z = Xpos, Ypos, Zpos 
PRINTER IS 16 
GOSUB Atod 


! CALCULATE SAMPLE VELOCITIES 

Ans$="N" 

INPUT "DO YOU WISH TO OBTAIN ESTIMATES OF U AND V ? CY/N, DEFAULT N)",Rns* 

IF Ans$="N" THEN GOTO 1720 

RAD 

ICALL F i nd_oe 1 
GOTO 2420 


WRITE CALIBRATION CONSTANTS AND DATA TO DISK 


GOSUB Df i 1 e 
Ans$="Y" 

INPUT "DO YOU WISH TO 
IF Ans$="N" THEN GOTO 
Run=Run+l 
Ans$="N" 

INPUT "DO YOU WISH TO 
IF Ans*="Y" THEN GOTO 


TAKE ANOTHER POINT ? CY/N, DEFAULT Y)",Ans$ 
1820 


CHANGE ANY PARAMETERS ? CY/N, DEFAULT N)",Ans$ 
540 
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1310 GOTO 1410 
1820 END 

1330 ! | m HM!!!!!!!!!!!!!!!!! M IM m mu m mm m | MMUMMM | 

1340 ! *********** END OF MAIN PROGRAM LDV ********************* 

1350 ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

I860 Atod: ! Subroutine for input from the LDV-A^D Cl 

1870 ! Enter routine with Ns = no. samples 

1880 DISP "Press CONT to initiate data acquisition" 

1390 PAUSE 

1900 DISP “Acquiring Data" 

1910 RESET 10 

1920 CONTROL MASK 10:; 1 

1930 WRITE 10 10,5;0 

1940 WRITE 10 10, 5; 1 ! start handshake by setting CTL0 

1950 Nt =6*Ns 

I960 FOR 1=1 TO 5 

1970 Dummy=READBIN<10> 

1980 NEXT I 

1990 REDIM Dl<Ns,6> 

2000 ! 

2010 ENTER 10 WFNS Nt NOFORMAT; D1 <*) Ifast data acquisition 

2020 WRITE 10 10,5;0 

2030 PRINT 

2048 DISP "Data acquisition complete" 

2050 PRINT 

2060 I CALL Data_tran:i 
2070 RETURN 

2080 ! !!!!) j !!!!! i !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! j !!!!!!!! j j !!!!!! ! 

2090 Dfile! ! write data file to floppy disc 

2100 PRINTER IS 16 
2110 PRINT 

2120 PRINT "** DATA FILE WRITE TO FLOPPY DISK **" 

2130 PRINT 

2140 PRINT "At this point be sure there is a floppy in drive 0 of" 

2150 PRINT "the 9895A with space for a file of 101, 256-byte records. " 

2160 PRINT 
2170 PRINT 
2180 Ans#="Y" 

2190 INPUT "DO YOU WISH TO WRITE THESE DATA TO DISK ? (Y/N, DEFAULT Y>",Ans# 
2200 IF Ans#="N" THEN GOTO 2380 
2210 Fi 1 e#=Name#&VAL# (Run) 

2220 DISP "File ";File#;" being written to disk" 

2230 MASS STORAGE IS ":H8,0,0" ! set floppy drive <9395R drive 0) as default 

2240 CREATE File#, 10:1 ! open file with 101 records 256 bytes each 

2250 ASSIGN Fi 1 e# TO #1 
2260 PRINT #1 ; Date# 

2270 PRINT # 1 ; T i t 1 # 

2280 PRINT #1 ; Name# 

2290 PRINT #1 ; Dn 

2300 PRINT #1 ; Nub, Numi xl , Numi x2, Numi x3, Theta, Run 
2310 PRINT #l;Vref,Ue,Dfl,Df2,Df3 
2320 PRINT #l;Xpos,Ypos,Zpos,Ns 
2330 MAT PRINT #l;Data 

2340 PRINT "***** File write completed *****" 

2350 ASSIGN * TO #1 ! close data file 

2360 MASS STORAGE IS ":H8,0,1" ! reset program disk as mass storage 

2370 GOTO 2390 
2380 Run=Run-l 
2390 RETURN 

2400 ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!( 
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2410 Uguess: ! Estimate U-component of velocity 

2420 DISP "CALCULATING ESTIMATES OF U AND V VELOCITIES" 

2430 ! 

2440 ! 

2450 IF Kill _w=3 THEN GOTO 2490 

2460 DATA 0,0,0 

2470 RESTORE 2460 

2480 READ Sw,Sww,Suw 

2490 Ubar=Su/Ns 

2500 Vbar=Sv/Ns 

2510 Wbar=SwxNs 

2520 Upr i 2=Suu/'Ns-Ubar*Ubar 

2530 Vpr i 2=Svv/'Ns-Vbar*Vbar 

2540 Wpr i 2=Sww/'Ns-Wbar*Wbar 

2550 Uvbar=Suv/'Ns-Ubar*Vbar 

2560 Uwbar=Suw/'Ns-Ubar*Wbar 

2570 Vwbar=Svw''Ns-Vbar*Wbar 

2580 PRINTER IS 0 

2590 PRINT "ESTIMATE OF QUANTITIES FROM SAMPLES : " , Ns 

2600 PRINT " Ubar =",Ubar 

2610 PRINT " Vbar =" ,Vbar 

2620 PRINT " Wbar =" , Nbar 

2630 PRINT "Upr i 2 =",Upri2 

2640 PRINT " Vpr i 2 =",Vpri2 

2650 PRINT "Wpr i 2 =",Wpri2 

2660 PRINT "Uvbar =", Uvbar 

2670 PRINT "Uwbar =", Uwbar 

2680 PRINT " Vwbar =", Vwbar 

2690 GOTO 1700 

2700 ISOURCE NAM Find_vel 

2710 ISOURCE ! 

2720 ISOURCE ! This subroutine converts raw data counts into instantaneous 
2730 ISOURCE ! velocities, then sums several different products of the 
2740 ISOURCE ! velocity components. All input and output data is passed 
2750 ISOURCE ! through the COMMON storage area. The inputs are the raw 
2760 ISOURCE ! data array <Arrayd), the Bragg shift frequency <Nub>, the 
2770 ISOURCE ! mixing frequencies CNumixl, Numix2, Numix3>, the fringe 
2780 ISOURCE ! spacings <Dfl, Df2, Df3), the crossing angle of the third 
2790 ISOURCE ! beam CTheta), and the number of samples in a data point <Ns>. 
2800 ISOURCE ! The outputs are the summations of various products of the 
2810 ISOURCE ! velocity components, including U, V, M, U*U, V*V, W*W, U*V, 
2820 ISOURCE ! U*W, and V*W. 

2830 ISOURCE ! 

2840 ISOURCE EXT Get_val ue ! Declare subroutines stored 

2850 ISOURCE EXT Get_info ! outside of the main program. 

2860 ISOURCE EXT Get_element 

2870 ISOURCE EXT Int_to_rel 

2880 ISOURCE EXT Re1_math 

2890 ISOURCE EXT Put_value 

2900 ISOURCE ! 

2910 ISOURCE COM 

2920 ISOURCE Data_par: INT <*> ! Declare common variables. 

2930 ISOURCE Dl_par: INT <*) 

2940 ISOURCE Ns_par: INT 

2950 ISOURCE Nn_par: INT 

2960 ISOURCE Dfl_par! REL 

2970 ISOURCE Df2_parl REL 

2980 ISOURCE Df3_parl REL 

2990 ISOURCE Theta_par: REL 

3000 ISOURCE Nub_par : REL 
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3010 

ISOURCE 

Nmi xl_par: 

REL 




3020 

ISOURCE 

Nm i x2_par : 

REL 




3030 

ISOURCE 

Nmi x3_par : 

REL 




3040 

ISOURCE 

Su_par : 

REL 




3050 

ISOURCE 

Su_par : 

REL 




3060 

ISOURCE 

Su_par : 

REL 




3070 

ISOURCE 

Suu_par : 

REL 




3080 

ISOURCE 

Suu_par : 

REL 




3090 

ISOURCE 

Suu_par :: 

REL 




3100 

ISOURCE 

Suu_par : 

REL 




3110 

ISOURCE 

Suu_par 3 

REL 




3120 

ISOURCE 

Suu_par l 

REL 




3130 

ISOURCE 

! 





3140 

ISOURCE 

Arrayd : 

BSS 

39 


! Reserue space for data array 

3150 

ISOURCE 

E 1 ementd : 

EQU 

Arrayd+16 ! descriptor. 

3160 

ISOURCE 

Array 1 : 

BSS 

4096 

! Reserue space for lookup tables used 

3170 

ISOURCE 

Array2: 

BSS 

4096 

! for count to uelocity conuersion. 

3180 

ISOURCE 

Array3: 

BSS 

4096 


3190 

ISOURCE 

Ns: 

BSS 

1 


! Reserue space for uarious input 

3200 

ISOURCE 

Df 1 : 

BSS 

4 


! and output uariables. 

3210 

ISOURCE 

Df2: 

BSS 

4 



3220 

ISOURCE 

Df3: 

BSS 

4 



3230 

ISOURCE 

Theta: 

BSS 

4 



3240 

ISOURCE 

Nub : 

BSS 

4 



3250 

ISOURCE 

Numi xl : 

BSS 

4 



3260 

ISOURCE 

Numi x2: 

BSS 

4 



3270 

ISOURCE 

Numi x3: 

BSS 

4 



3280 

ISOURCE 

Su: 

BSS 

4 



3290 

ISOURCE 

Su: 

BSS 

4 



3300 

ISOURCE 

Su: 

BSS 

4 



3310 

ISOURCE 

Suu: 

BSS 

4 



3320 

ISOURCE 

Suu: 

BSS 

4 



3330 

ISOURCE 

Suu: 

BSS 

4 



3340 

ISOURCE 

Suu: 

BSS 

4 



3350 

ISOURCE 

Suu: 

BSS 

4 



3360 

ISOURCE 

Suu: 

BSS 

4 



3370 

ISOURCE 

Count : 

BSS 

1 

| 

Count and I are general purpose index uari- 

3380 

ISOURCE 

i: 

BSS 

1 

! 

ables. Count is usually 0, 1, or 2, to denote 

3390 

ISOURCE 

Check : 

BSS 

1 

i 

uhether U,V, or W is being calculated. 

3400 

ISOURCE 

Int : 

BSS 

1 

1 

Int, Address, and Offset are all general 

3410 

ISOURCE 

Address;! 

BSS 

1 

i 

purpose storage areas. 

3420 

ISOURCE 

Offset : 

BSS 

1 



3430 

ISOURCE 

Rl: 

BSS 

4 

j 

Rl, R2, and R3 are the count-to- 

3440 

ISOURCE 

R2: 

BSS 

4 

I 

frequency conuersion factors. 

3450 

ISOURCE 

R3: 

BSS 

4 



3460 

ISOURCE 

Xuar : 

BSS 

4 

i 

Xuar and Yuar are general purpose real 

3470 

ISOURCE' 

Yuar: 

BSS 

4 

j 

number storage areas. 

3480 

ISOURCE 

u: 

BSS 

4 

! 

reserue space for instantaneous uelocity 

3490 

ISOURCE 

V: 

BSS 

4 

j 

c omponent s . 

3500 

ISOURCE 

w: 

BSS 

4 



3510 

ISOURCE 

Uu: 

BSS 

4 



3520 

ISOURCE 

uu: 

BSS 

4 



3530 

ISOURCE 

Cos: 

BSS 

4 

! 

Cos and Sin are the cosine and sine of Theta. 

3540 

ISOURCE 

Sin: 

BSS 

4 



3550 

ISOURCE 

Rad: 

DAT 

5. 

729578E 1 

3560 

ISOURCE 

Mill: 

DAT 

1. 

E6 


3570 

ISOURCE 

One: 

DAT 

1. 



3580 

ISOURCE 

Zero: 

DAT 

0. 



3590 

ISOURCE 

! 





3600 

ISOURCE 


SUB 
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3610 

ISOURCE 

F i nd_ 

ye 1 : 

LDfl 

=Ns ! 

Get number of samples. 

3620 

ISOURCE 



LDB 

=Ns_par 


3630 

ISOURCE 



JSM 

Get_ya1 ue 


3640 

ISOURCE 



LDfl 

=Arrayd ! 

Get parameters of data array. 

3650 

ISOURCE 



LDB 

=Data_par 


3660 

ISOURCE 



JSM 

Get info 


3670 

ISOURCE 



LDfl 

=Df 1 ! 

Get input parameters. 

3680 

ISOURCE 



LDB 

=Df l_par 


3690 

ISOURCE 



JSM 

Get yal ue 


3700 

ISOURCE 



LDfl 

=Df 2 


3710 

ISOURCE 



LDB 

=Df2_par 


3720 

ISOURCE 



JSM 

Get yal ue 


3730 

ISOURCE 



LDfl 

= Df 3 


3740 

ISOURCE 



LDB 

=Df 3_par 


3750 

ISOURCE 



JSM 

Get_ual ue 


3760 

ISOURCE 



LDfl 

=Theta 


3770 

ISOURCE 



LDB 

=Thet a_par 


3780 

ISOURCE 



JSM 

Get_yal ue 


3790 

ISOURCE 



LDfl 

=Nub 


3800 

ISOURCE 



LDB 

=Nub_par 


3810 

ISOURCE 



JSM 

Get_val ue 


3820 

ISOURCE 



LDfl 

=Numi xl 


3830 

ISOURCE 



LDB 

=Nmi xl_par 


3840 

ISOURCE 



JSM 

Get_va1 ue 


3850 

ISOURCE 



LDfl 

=Numi x2 


3860 

ISOURCE 



LDB 

=Nm i x2_par 


3870 

ISOURCE 



JSM 

Get_ya1 ue 


3880 

ISOURCE 



LDfl 

=Numi x3 


3890 

ISOURCE 



LDB 

=Nmi x3_par 


3900 

ISOURCE 



JSM 

Get_val ue 


3910 

ISOURCE 

! 





3920 

ISOURCE 

! The 

1 oop 

headed by Get free 

is repeated three times to get 

3930 

ISOURCE 

! the 

count 

-to- 

■frequency conversion factors (which depend on 

3940 

ISOURCE 

! the 

range) for U, V, and W. 

Whenever a loop is controlled by 

3950 

ISOURCE 

! the 

vari abl e 

"Count", the 1 

oop contains operations which are 

3960 

ISOURCE 

! the 

same 

for 

U, V, and W. 


3970 

ISOURCE 

j 





3980 

ISOURCE 



LDfl 

= 0 


3990 

ISOURCE 



STfl 

Count 


4000 

ISOURCE 

Get freq: 

LDfl 

Count ! 

Get the first word of the column of 

4010 

ISOURCE 



LDB 

Ns ! 

the data array which contains the 

4020 

ISOURCE 



MPY 

1 

velocity component for which we want 

4030 

ISOURCE 



STfl 

Elementd ! 

to get the range. 

4040 

ISOURCE 



LDfl 

= Int 


4050 

ISOURCE 



LDB 

=firrayd 


4060 

ISOURCE 



JSM 

Get element 


4070 

ISOURCE 



LDfl 

Int ! 

Mask and rotate to get the four 

4080 

ISOURCE 



LDB 

=15360 ! 

bits containing the range. 

4090 

ISOURCE 



AND 

B 


4100 

ISOURCE 



SflR 

10 


4110 

ISOURCE 



TCfl 

; 

Subtract from 15 to get the 

4120 

ISOURCE 



LDB 

= 15 ! 

actual range. 

4130 

ISOURCE 



ADA 

B 


4140 

ISOURCE 



LDB 

= 1 


4150 

ISOURCE 



SZfl 

Loopend 


4160 

ISOURCE 

Loop: 


SBL 

1 ! 

Use the range to find the power 

4170 

ISOURCE 



DSZ 

fl ! 

of two needed for the divisor. 

4180 

ISOURCE 



JMP 

Loop 


4190 

ISOURCE 

Loopend : 

STB 

Int 


4200 

ISOURCE 



LDfl 

= I nt ! 

Convert the power of two into a 
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4210 

ISOURCE 

STfl 

Oper_l 

4220 

ISOURCE 

LDfl 

=Yvar 

4230 

ISOURCE 

STfl 

Resul t 

4240 

ISOURCE 

JSM 

Int to rel 

4250 

ISOURCE 

LDfl 

==3 . 2E4 

4260 

ISOURCE 

LDB 

=Xvar 

4270 

ISOURCE 

XFR 

4 

4280 

ISOURCE 

STB 

Oper_l 

4290 

ISOURCE 

LDfl 

=Yvar 

4300 

ISOURCE 

STfl 

0per_2 

4310 

ISOURCE 

LDfl 

Count 

4320 

ISOURCE 

LDB 

=4 

4330 

ISOURCE 

MPY 


4340 

ISOURCE 

ADA 

=R 1 

4350 

ISOURCE 

STfl 

Resul t 

4360 

ISOURCE 

LDfl 

= 2 

4370 

ISOURCE 

LDB 

■ 1 47 1 55B 

4380 

ISOURCE 

JSM 

Rel_math 

4390 

ISOURCE 

ISZ 

Count 

4480 

ISOURCE 

LDfl 

= 3 

4410 

ISOURCE 

CPA 

Count 

4420 

ISOURCE 

JMP 

*+2 

4430 

ISOURCE 

JMP 

Get_f req 

4440 

ISOURCE ! 



4450 

ISOURCE 

LDfl 

=Array 1 

4460 

ISOURCE 

LDB 

=768 

4470 

ISOURCE Continue: 

CLR 

16 

4480 

ISOURCE- 

ADA 

= 16 

4490 

ISOURCE 

DSZ 

B 

4500 

ISOURCE 

JMP 

Cont i nue 

4510 

ISOURCE ! 



4520 

ISOURCE 

LDfl 

=Su 

4530 

ISOURCE 

LDB 

=9 

4540 

ISOURCE Clear: 

CLR 

4 

4550 

ISOURCE 

ADA 

=4 

4560 

ISOURCE 

DSZ 

B 

4570 

ISOURCE 

JMP 

Clear 

4580 

ISOURCE ! 



4590 

ISOURCE 

LDfl 

= Thet a 

4600 

ISOURCE 

STfl 

Oper 1 

4610 

ISOURCE 

LDfl 

= Rad 

4620 

ISOURCE 

STfl 

0per_2 

4630 

ISOURCE 

LDfl 

=Xvar 

4640 

ISOURCE 

STfl 

Resul t 

4650 

ISOURCE 

LDfl 

=2 

4660 

ISOURCE 

LDB 

=147155B 

4670 

ISOURCE 

JSM 

Rel_m ath 

4680 

ISOURCE ! 



4690 

ISOURCE 

LDfl 

=Xvar 

4700 

ISOURCE 

STfl 

Oper_l 

4710 

ISOURCE 

LDfl 

=S i n 

4720 

ISOURCE 

STfl 

Resul t 

4730 

ISOURCE 

LDfl 

= 1 

4740 

ISOURCE 

LDB 

=342 13B 

4750 

ISOURCE 

JSM 

Rel_math 

4760 

ISOURCE 

LDfl 

= Cos 

4770 

ISOURCE 

STfl 

Resul t 

4780 

ISOURCE 

LDfl 

= 1 

4790 

ISOURCE 

LDB 

=34224B 

4800 

ISOURCE , 

JSM 

Re l_math 


! real number. 


! Divide 3.2E4 by the appropriate 
! power of two, using BCD math. 


Decide whether to put the result 
in Rl, R2, or R3, depending on 
Count . 


! Now, finally, call the utility to 
! perform the division. 

! Increment and check Count so as 
! to follow the loop three times. 


! Zero out the entire count-to- 
! velocity conversion table so that 
! it must be recalculated for each 
! point. (This must be done if the 
! mixing frequencies or ranges are 
! changed between counts. > 

! Set initial values of Su, Sv, 

! Sw, Suu, etc. to zero. 


Convert Theta from degrees to 
radians using the Rel_math 
ut i 1 i ty . 


Find the sine and cosine of 
Theta, and store them in the 
locations Sin and Cos, respec t i ve 1 y . 
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4810 

ISOURCE 

! 



4820 

ISOURCE 

l 

The loop defined by Get_int . cal culates intermediate values 

4830 

ISOURCE 

j 

used in converting counts to velocities. These values are the 

4840 

ISOURCE 

j 

same for al 1 

samples in a point, so they can be calculated 

4350 

ISOURCE 

j 

separately. The loop calculates (Nub-NumixN) and <Mi11*DfN), 

4860 

ISOURCE 

j 

where N is 1, 

2, and 3. 

4870 

ISOURCE 

j 



4880 

ISOURCE 


LDfl 

=3 

4890 

ISOURCE 


STfl 

Count 

4900 

ISOURCE 

Get int:: LDfl 

Count 

4910 

ISOURCE 


ADA 

=-l 

4920 

ISOURCE 


SflL 

2 

4930 

ISOURCE 


STfl 

Offset 

4940 

ISOURCE 


LDfl 

=Nub ! Find Num i xN=Nub-Num i xN . 

4950 

ISOURCE 


STfl 

Oper_l 

4960 

ISOURCE 


LDfl 

=Numi xl 

4970 

ISOURCE 


ADA 

Offset 

4980 

ISOURCE 


STfl 

0per_2 

4990 

ISOURCE 


STfl 

Resul t 

5000 

ISOURCE 


LDfl 

= 2 

5010 

ISOURCE 


LDB 

=146717B 

5020 

ISOURCE 


JSM 

Rel math 

5030 

ISOURCE 


LDfl 

=Df 1 ! Find Df N=Mi 1 1 *Df N. 

5040 

ISOURCE 


ADA 

Offset 

5050 

ISOURCE 


STfl 

Oper_l 

5060 

ISOURCE 


STfl 

Result 

5070 

ISOURCE 


LDfl 

=Mi 1 1 

5080 

ISOURCE 


STfl 

Oper 2 

5090 

ISOURCE 


LDfl 

=2 

5100 

ISOURCE 


LDB 

= 1 47037B 

5110 

ISOURCE 


JSM 

Rel_math 

5120 

ISOURCE 


DSZ 

Count 

5130 

ISOURCE 


JMP 

Get_i nt 

5140 

ISOURCE 

J 



5150 

ISOURCE 

J 

The loop Begin is performed three times, once for U, V, 

5160 

ISOURCE 

1 

and W. Each time through the data array is read, a count <a 

5170 

ISOURCE 

| 

raw datum) is 

taken from it, and is converted into a velocity. 

5180 

ISOURCE 

J 

The velocity 

is stored in U, V, or W depending on whether 

5190 

ISOURCE 

j 

this is the first, second, or third iteration of the loop. The 

5200 

ISOURCE 

j 

first time any particular count is encountered, the velocity 

5210 

ISOURCE 

J 

correspondi ng to it is calculated using the intermediate values 

5220 

ISOURCE 

j 

found i n Get 

int, and the velocity is stored in a table. If that 

5230 

ISOURCE 

I 

count is found again in the data array, the c orrespond i ng vel- 

5240 

ISOURCE 

J 

ocity is looked up rather than being calculated again. 

5250 

ISOURCE 

J 



5260 

ISOURCE 


LDfl 

=0 

5270 

ISOURCE 


STfl 

Count 

5280 

ISOURCE 


LDfl 

=5 

5290 

ISOURCE 


STfl 

Check 

5380 

ISOURCE 


LDfl 

Ns 

5310 

ISOURCE 


STfl 

I 

5320 

ISOURCE 

! 



5330 

ISOURCE 

Begin: LDfl 

Count 

5340 

ISOURCE 


STfl 

B 

5350 

ISOURCE 


SBL 

2 

5360 

ISOURCE 


STB 

Offset 

5370 

ISOURCE 


LDB 

Ns ! Figure out which element of the 

5380 

ISOURCE 


MPY 

! data array we want to pick up. 

5390 

ISOURCE 


ADA 

I 

5400 

ISOURCE 


ADA 

=-l 
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5410 

ISOURCE 


STfl 

El ementd 

5420 

ISOURCE 


LDfl 

= I nt 

5430 

ISOURCE 


LDB 

=flrrayd 

5440 

ISOURCE 


JSM 

Get_el ement 

5450 

ISOURCE 


LDfl 

Int 

5460 

ISOURCE 


LDB 

= 1023 

5470 

ISOURCE 


FIND 

B 

5480 

ISOURCE 


LDB 

Count 

5490 

ISOURCE 


TCB 


5500 

ISOURCE 


flDB 

Check 

5510 

ISOURCE 


SZB 

Strai ght 

5520 

ISOURCE 


TCfi 


5530 

ISOURCE 


ADA 

= 1024 

5540 

ISOURCE 

St rai ght : 

STfl 

Int 

5550 

ISOURCE 


LDfl 

Count 

5560 

ISOURCE 


LDB 

=4096 

5570 

ISOURCE 


MPY 


5580 

ISOURCE 


LDB 

Int 

5590 

ISOURCE 


flDB 

= -l 

5600 

ISOURCE 


SBL 

2 

5610 

ISOURCE 


ADA 

B 

5620 

ISOURCE 


ADA 

= 1 

5630 

ISOURCE 


ADA 

=flrray 1 

5640 

ISOURCE 


STfl 

Address 

5650 

ISOURCE 


LDfl 

Address, I 

5660 

ISOURCE 


SZfl 

Cal cul ate 

5670 

ISOURCE 


JMP 

Over 

5680 

ISOURCE 

Cal cul a>e: 

LDfl 

= I nt 

5690 

ISOURCE 


STfl 

Oper_l 

5700 

ISOURCE 


LDfl 

= Yvar 

5710 

ISOURCE 


STfl 

Resul t 

5720 

ISOURCE 


JSM 

Int_t o_rel 

5730 

ISOURCE 


STB 

Oper 2 

5740 

ISOURCE 


LDfl 

= R 1 

5750 

ISOURCE 


ADA 

Offset 

5760 

ISOURCE 


STfl 

Oper_l 

5770 

ISOURCE 


LDfl 

=Xvar 

5780 

ISOURCE 


STfl 

Resul t 

5790 

ISOURCE 


LDfl 

=2 

5800 

ISOURCE 


LDB 

=1471 55B 

5810 

ISOURCE 


JSM 

Rel_math 

5820 

ISOURCE 

1 



5830 

ISOURCE 


LDfl 

=Num i x 1 

5840 

ISOURCE 


ADA 

Offset 

5850 

ISOURCE 


STfl 

Oper_l 

5860 

ISOURCE 


LDfl 

=Xvar 

5870 

ISOURCE 


STfl 

0per_2 

5880 

ISOURCE 


LDfl 

=Yvar 

5890 

ISOURCE 


STfl 

Resul t 

5900 

ISOURCE 


LDfl 

=2 

5910 

ISOURCE 


LDB 

=146717B 

5920 

ISOURCE 


JSM 

Re l_mat h 

5930 

ISOURCE 

! 



5940 

ISOURCE 


LDB 

Count 

5950 

ISOURCE 


flDB 

=-l 

5960 

ISOURCE 


SBM 

*+2 

5970 

ISOURCE 


JMP 

Samesi gn 

5980 

ISOURCE 


LDfl 

= Zero 

5990 

ISOURCE 


STfl 

Oper_l 

6000 

ISOURCE 


LDfl 

Resul t 


! Get a raw datum from the data array. 


Strip off the first six 
bits of the raw data word. 
See if Count = Check. 


If true, use the stripped data 
word as an index. If not, use 
1024 minus the data word. 

Store the count we have gotten. 
Now, use Count to find out 
which lookup table array 
we want to use, and use the 
count we got from the data array 
to find exactly where in the table 
we want to go. 


! If that table entry is zero, 
! calculate a velocity for it. 


! Convert the count to a real number. 


Divide the range we found 
earlier by the count to get a 
f requency . 


! Find <Nub-Numi xN)-FrequencyN. 


If we are calculating U, reverse 
the sign of < <Nub-Numi xN)-FrequencyhD 
so as to reverse the sign of U. 

Leave V and M alone. 
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6010 

ISOURCE 

STfl 

0per_2 

6020 

ISOURCE 

LDfl 

=Xvar 

6030 

ISOURCE 

STfl 

Resu 1 t 

6040 

ISOURCE 

LDfl 

= 2 

6050 

ISOURCE 

LDB 

= 1 467 1 7B 

6060 

ISOURCE 

JSM 

Rel_math 

6070 

ISOURCE 

Samesign: NOP 


6080 

ISOURCE 

! 


6090 

ISOURCE 

LDfl 

Resul t 

6100 

ISOURCE 

STfl 

Oper 1 

6110 

ISOURCE 

LDfl 

=Df 1 

6120 

ISOURCE 

ADA 

Offset 

6130 

ISOURCE 

STfl 

0per_2 

6140 

ISOURCE 

LDfl 

Address 

6150 

ISOURCE 

ADA 

=-l 

6160 

ISOURCE 

STfl 

Resul t 

6170 

ISOURCE 

LDfl 

=2 

6180 

ISOURCE 

LDB 

= 1 47037B 

6190 

ISOURCE 

JSM 

Rel_math 

6200 

ISOURCE 

! 


6210 

ISOURCE 

Over: LDfl 

Address 

6220 

ISOURCE 

ADA 

= - 1 

6230 

ISOURCE 

LDB 

=U 

6240 

ISOURCE 

ADB 

Offset 

6250 

ISOURCE 

XFR 

4 

6260 

ISOURCE 

! 


6270 

ISOURCE 

ISZ 

Count 

6280 

ISOURCE 

LDfl 

=-3 

6290 

ISOURCE 

ADA 

Count 

6300 

ISOURCE 

Pause: HOP 


6310 

ISOURCE 

SZfl 

*+2 

6320 

ISOURCE 

JMP 

Begi n 

6330 

ISOURCE 

STfl 

Count 

6340 

ISOURCE 

! 


6350 

ISOURCE 

! Now we convert the w w 

6360 

ISOURCE 

! an angle Theta to the 

6370 

ISOURCE 

! be measured 

at an angl 

6380 

ISOURCE 

! Thus find W= 

(Wv-V*Cos( 

6390 

ISOURCE 

! 


6400 

ISOURCE 

LDfl 

=V 

6410 

ISOURCE 

STfl 

Oper_l 

6420 

ISOURCE 

LDfl 

=Cos 

6430 

ISOURCE 

STfl 

0per_2 

6440 

ISOURCE 

LDfl 

=Xvar 

6450 

ISOURCE 

STfl 

Resul t 

6460 

ISOURCE 

LDfl 

= 2 

6470 

ISOURCE 

LDB 

= 1 47037B 

6480 

ISOURCE 

JSM 

Rel_math 

6490 

ISOURCE 

LDfl 

Resu 1 t 

6500 

ISOURCE 

STfl 

Oper 2 

6510 

ISOURCE 

LDfl 

=W 

6520 

ISOURCE 

STfl 

Oper_l 

6530 

ISOURCE 

LDfl 

=Yvar 

6540 

ISOURCE 

STfl 

Resul t 

6550 

ISOURCE 

LDfl 

=2 

6560 

ISOURCE 

LDB 

= 1 467 1 7B 

6570 

ISOURCE 

JSM 

Rel_math 

6580 

ISOURCE 

LDfl 

Resul t 

6590 

ISOURCE 

STfl 

Oper_l 

6600 

ISOURCE 

LDfl 

= Si n 


Find Vel oc i ty=< (Nub-Numi xN) 
-FrequencyN;>*(Mi 1 1 *DfN> 
and store in a place in the 
lookup table correspondi ng to 
the data count. 


How transfer the velocity from 
the lookup table to U, V, or W, 
as appropriate. 


! Have U, V, and W al 1 
! been calculated? 

! If not, go back again. 

! If so, set Count = 8. 

e have obtained (which is measured at 
V-axis) to the W we want (which should 
e of 90 degrees to the V-axis>. 
Theta))xSin(Theta!J. 
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6610 

ISOURCE 



STfl 

Oper 2 



6620 

ISOURCE 



LDfl 

= W 



6630 

ISOURCE 



STfl 

Result 



6640 

ISOURCE 



LDfl 

=2 



6650 

ISOURCE 



LOB 

= 1 47 1 55B 



6660 

ISOURCE 



JSM 

Rel_math 



6670 

ISOURCE 

! 






6680 

ISOURCE 

j 

Now take 

running sums of U 

, V, W, and several 

produc t s 

6690 

ISOURCE 

! 

of these 

velocities. The sums are taken using 

the utility 

6700 

ISOURCE 

! 

"Add", The sums are calculated in an unusual 

sequence in 

6710 

ISOURCE 

! 

order to 

reduce the number 

of steps needed to 

cal cul ate 

6720 

ISOURCE 

j 

them. Program steps which 

produce sums not considered 

6730 

ISOURCE 

! 

necessary have been included, but these steps 

are preceded 

6740 

ISOURCE 

! 

by the "< 

" comment marker. 



6750 

ISOURCE 

! 






6760 

ISOURCE 



LDfl 

= Su ! 

Find Su=Su+U 


6770 

ISOURCE 



LOB 

=U 



6780 

ISOURCE 



JSM 

Add 



6790 

ISOURCE 

! 






6880 

ISOURCE 



LDfl 

0per_2 ! 

Find U*U 


6810 

ISOURCE 



STfl 

Oper 1 



6820 

ISOURCE 



LDfl 

=Uu 



6830 

ISOURCE 



STfl 

Resul t 



6840 

ISOURCE 



LDfl 

= 2 



6850 

ISOURCE 



LDB 

=147037B 



6860 

ISOURCE 



JSM 

Rel_math 



6870 

ISOURCE 

! 






6880 

ISOURCE 



LDfl 

=Suu ! 

Find Su=Su+<U*U> 


6890 

ISOURCE 



LDB 

Resul t 



6900 

ISOURCE 



JSM 

Add 



6910 

ISOURCE 

! 






6920 

ISOURCE 



LDfl 

=Su ! 

Find Sv=Sv+V 


6930 

ISOURCE 



LDB 

=V 



6940 

ISOURCE 



JSM 

Add 



6950 

ISOURCE 

! 






6960 

ISOURCE 



LDfl 

=U ! 

Find U*V 


6970 

ISOURCE 



STfl 

Oper 1 



6980 

ISOURCE 



LDfl 

=Uv 



6990 

ISOURCE 



STfl 

Resu 1 1 



7000 

ISOURCE 



LDfl 

=2 



7010 

ISOURCE 



LDB 

=147037B 



7020 

ISOURCE 



JSM 

Rel_math 



7030 

ISOURCE 

! 






7040 

ISOURCE 



LDfl 

=Suv ! 

Find Suv=Suv+ < U*V > 


7050 

ISOURCE 



LDB 

Resul t 



7060 

ISOURCE 



JSM 

Add 



7070 

ISOURCE 

! 






7080 

ISOURCE 



LDfl 

=V ! 

Find 


7090 

ISOURCE 



STfl 

Oper_l 



7100 

ISOURCE 



STfl 

0per_2 



7110 

ISOURCE 



LDfl 

=Xvar 



7120 

ISOURCE 



STfl 

Resul t 



7130 

ISOURCE 



LDfl 

=2 



7140 

ISOURCE 



LDB 

= 147037B 



7150 

ISOURCE 



JSM 

Re 1 _mat h 



7160 

ISOURCE 

! 






7170 

ISOURCE 



LDfl 

=Svv ! 

Find Svv=Svv+< V*V) 


7180 

ISOURCE 



LDB 

Resul t 



7190 

ISOURCE 



JSM 

Add 



7200 

ISOURCE 

! 
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7210 

ISOURCE 



LDfl 

= Sw 

! Find 

S w = S w + M 

7220 

ISOURCE 



LDB 

= kl 



7230 

ISOURCE 



JSM 

Add 



7240 

ISOURCE 

J 






7250 

ISOURCE 



LDfl 

= U 

! Find 

U*W 

7260 

ISOURCE 



STfl 

Oper_l 



7270 

ISOURCE 



LDfl 

=Xvar 



7280 

ISOURCE 



STfl 

Resul t 



7290 

ISOURCE 



LDfl 

=2 



7300 

ISOURCE 



LDB 

=147037B 



7310 

ISOURCE 



JSM 

Rel_math 



7320 

ISOURCE 

! 






7330 

ISOURCE 



LDfl 

=Suw 

! Find 

Suw=Suw+(U*W> 

7340 

ISOURCE 



LDB 

Resul t 



7350 

ISOURCE 



JSM 

Add 



7360 

ISOURCE 

! 






7370 

ISOURCE 



LDfl 

=W 

! Find 

V*W 

7380 

ISOURCE 



STfl 

Oper 1 



7390 

ISOURCE 



LDfl 

=V 



7400 

ISOURCE 



STfl 

0per_2 



7410 

ISOURCE 



LDfl 

=Xvar 



7420 

ISOURCE 



STfl 

Resul t 



7430 

ISOURCE 



LDfl 

= 2 



7440 

ISOURCE 



LDB 

= 1 47037B 



7450 

ISOURCE 



JSM 

Rel_math 



7460 

ISOURCE 

! 






7470 

ISOURCE 



LDfl 

=Svw 

! Find 

Svw=Svw+< V*W) 

7480 

ISOURCE 



LDB 

Resul t 



7490 

ISOURCE 



JSM 

Add 



7500 

ISOURCE 

! 






7510 

ISOURCE 



LDfl 

=W 

! Find 

w*w 

7520 

ISOURCE 



STfl 

Oper_l 



7530 

ISOURCE 



STfl 

Gper_2 



7540 

ISOURCE 



LDfl 

=Xvar 



7550 

ISOURCE 



STfl 

Resul t 



7560 

ISOURCE 



LDfl 

=2 



7570 

ISOURCE 



LDB 

= 1 47037B 



7580 

ISOURCE 



JSM 

Rel_math 



7590 

ISOURCE 

j 






7600 

ISOURCE 



LDfl 

= Sww 

! Find 

Sww=Sww+<W*W> 

7610 

ISOURCE 



LDB 

Resul t 



7620 

ISOURCE 



JSM 

Add 



7630 

ISOURCE 

! 






7640 

ISOURCE 



DS2 

I 

! Continue to calculate running 

7650 

ISOURCE 



JMP 

Begi n 

! sums 

until out of samples. 

7660 

ISOURCE 

j 






7670 

ISOURCE 

j 

Place the finished sums 

in the 

COMMON region so that 

7680 

ISOURCE 

i 

the 

BASIC program has access to 

them, and then return to 

7690 

ISOURCE 

! 

the 

BASIC program. 



7700 

ISOURCE 

i 






7710 

ISOURCE 



LDfl 

=Su 



7720 

ISOURCE 



LDB 

=Su_par 



7730 

ISOURCE 



JSM 

Put value 



7740 

ISOURCE 



LDfl 

=Sv 



7750 

ISOURCE 



LDB 

=Sv_par 



7760 

ISOURCE 



JSM 

Put_val ue 



7770 

ISOURCE 



LDfl 

=Sw 



7730 

ISOURCE 



LDB 

=Sw_par 



7790 

ISOURCE 



JSM 

Put_val ue 



7800 

ISOURCE 



LDfl 

=Suu 
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7310 

ISOURCE 


LDB 

=Suu_par 



7820 

ISOURCE 


JSM 

Put_va1 ue 



7830 

ISOURCE 


LDA 

=Svv 



7840 

ISOURCE 


LDB 

=Svv_par 



7850 

ISOURCE 


JSM 

Put_oal ue 



7860 

ISOURCE 


LDA 

= Suu 



7870 

ISOURCE 


LDB 

=Sww_par 



7880 

ISOURCE 


JSM 

Put_va1 ue 



7890 

ISOURCE 


LDA 

=Suu 



7900 

ISOURCE 


LDB 

=Suv_par 



7910 

ISOURCE 


JSM 

Put_oa1 ue 



7920 

ISOURCE 


LDA 

=Suu 



7930 

ISOURCE 


LDB 

=Suu_par 



7940 

ISOURCE 


JSM 

Put_oal ue 



7950 

ISOURCE 


LDA 

=Svu 



7960 

ISOURCE 


LDB 

=Svw_par 



7970 

ISOURCE 


JSM 

Put value 



7980 

ISOURCE 


RET 

1 



7990 

ISOURCE 

! 





8000 

ISOURCE 

! The ut i 1 

i ty ' 

Add" is used 

to add up 

the running sums. 

8010 

ISOURCE 

! 





8020 

ISOURCE 

Add: 

IS2 

Ut 1 count 



8030 

ISOURCE 


STA 

Oper_l 



8040 

ISOURCE 


STB 

0per_2 



8050 

ISOURCE 


STA 

Resul t 



8060 

ISOURCE 


LDA 

=2 



8070 

ISOURCE 


LDB 

=146721B 



8080 

ISOURCE 


JSM 

Re1_math 



8090 

ISOURCE 


DSZ 

Ut 1 count 



8100 

ISOURCE 


RET 

1 



8110 

ISOURCE 


JSM 

Ut 1 end 



8120 

ISOURCE 

! 





8130 

ISOURCE 


LIT 

200 



8140 

ISOURCE 


END 

F i nd_vel 



8150 

ISOURCE 

! 





8160 

ISOURCE 


NAM 

Data_trans 



8170 

ISOURCE 

! 





8180 

ISOURCE 

! Subrout i ne tc 

transfer raw 

data from input array 

8190 

ISOURCE 

! to storage array. 



8200 

ISOURCE 

! 





8210 

ISOURCE 


EXT 

Get_info ! 

Declare 

subroutines stored 

8220 

ISOURCE 


EXT 

Get_value ! 

outs i de 

of the main program 

8230 

ISOURCE 


EXT 

Get_e1 ement 



8240 

ISOURCE 


EXT 

Put_e 1 ement 



8250 

ISOURCE 

i 





8260 

ISOURCE 


COM 




8270 

ISOURCE 

D at a_p ar : 

INT 

<*) ! 

Dec 1 are 

common variables. 

8280 

ISOURCE 

D 1 _p ar : 

INT 

<*> 



8290 

ISOURCE 

Ns_par : 

INT 




8300 

ISOURCE 

! 





8310 

ISOURCE 

Arrayd: 

BSS 

39 



8320 

ISOURCE 

El ementd: 

EQU 

Arrayd+ 1 6 



8330 

ISOURCE 

Array 1 : 

BSS 

39 



8340 

ISOURCE 

El ement 1 : 

EQU 

Array 1 + 16 



8350 

ISOURCE 

Ns: 

BSS 

1 



8360 

ISOURCE 

Count 1 : 

BSS 

1 



8370 

ISOURCE 

Count2: 

BSS 

1 



8380 

ISOURCE 

I: 

BSS 

1 



8390 

ISOURCE 

Int : 

BSS 

1 



8400 

ISOURCE 

! 
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8410 

8420 

8430 

8440 

8450 

8460 

8470 

8480 

8490 

8500 

8510 

8520 

8530 

8540 

8550 

8560 

8570 

8580 

8590 

8600 

8610 

8620 

8630 

8640 

8650 

8660 

8670 

8630 

8690 

8700 

8710 

8720 

8730 

8740 

8750 

8760 

8770 

8780 

8790 

8800 

8810 

8820 

8830 

8840 

8850 

8860 

8870 

8880 

8890 

8900 

8910 

8920 

8930 

8940 

8950 


ISOURCE 

I SOURCE Data_trans 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE ! 

ISOURCE 

ISOURCE 

ISOURCE St art 1 : 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE ! 

ISOURCE St art 2 : 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE ! 

ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE 
ISOURCE ! 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 

ISOURCE 


SUB 

LDA =Arrayd 


LDB 

=Data_par 

JSM 

Get_i nf o 

LDA 

=Array 1 

LDB 

=Dl_par 

JSM 

Get info 

LDA 

=Ns 

LDB 

=Ns_par 

JSM 

Get_oal ue 

LDA 

= 0 

STA 

I 

LDA 

= 0 

STA 

Count 1 

LDA 

=2 

STA 

Count2 

LDA 

I 

LDB 

=6 

MPY 

ADA 

Count2 

STA 

El ement 1 

LDA 

= Int 

LDB 

=Array 1 

JSM 

Get_el ement 

LDA 

Count 1 

LDB 

Ns 

MPY 

ADA 

I 

STA 

Elementd 

LDA 

= Int 

LDB 

=Rrrayd 

JSM 

Put_el ement 

ISZ 

Count 1 

ISZ 

Count2 

LDA 

=-3 

ADA 

Count 1 

SZA 

*+2 

JMP 

St art 2 

ISZ 

I 

LDA 

I 

TCA 

ADA 

Ns 

SZA 

*+2 

JMP 

Start 1 

RET 

1 

END 

Data trans 


Get parameters of data storage array. 
Get parameters of data input array. 
Get number of samples. 


************************************************************************ 
** HISTOGRAM SUBPROGRAM ** 

************************************************************************ 


8960 Hist: ! Subroutine to produce online 

8970 ! histograms of raw LBV data. 

8980 ON KEY #0 GOTO 9270 

8990 N n = 0 

9000 Ns_temp=Ns 
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9010 ! 

9020 INPUT "No. of data samples per point <2000 maxm. ) : " , Ns 
9030 ! 

9040 REDIM Data<3, Ns > 

9850 REDIM Dl<Ns,6) 

9060 DISP "Press CONT to initiate data acquisition, press K0 to return to main 
program. " 

9070 PAUSE 
9080 GCLEAR 
9090 GRAPHICS 
9100 ! 

9110 RESET 10 

9120 CONTROL MASK 105 1 

9130 WRITE 10 10, 5; 0 

9140 WRITE 10 10,5; 1 ! start handshake by settina CTL0 

9150 Nt=6*Ns 

9160 FOR 1=1 TO 5 

9178 Dummy=READBIN<10> 

9180 NEXT I 

9190 ! 

9200 ENTER 10 WFHS Nt NOFORMAT; D1 <*> !fast data acquisition 

9210 WRITE 10 1 0 . 5 ; 0 

9220 ICALL D4ta_trans 

9230 ICALL Draw_hist 

9240 IF Ns=2000 THEN GOTO 9270 

9250 Nn= 1 

9260 GOTO 9110 

9270 IF Ns=2000 THEN DUMP GRAPHICS 
9280 WRITE 10 10,5;0 
9290 Ns=Ns_temp 
9300 REDIM Data(3,Ns> 

9310 REDIM Dl(Ns,6> 

9320 EXIT GRAPHICS 
9330 RETURN 
9340 ! 

9350 ISOURCE NAM Draw_hist 

9360 ISOURCE 
9370 ISOURCE 
9380 ISOURCE 
9390 ISOURCE 
9400 ISOURCE 
9410 ISOURCE 
9420 ISOURCE 
9430 ISOURCE 
9440 ISOURCE 
9450 ISOURCE 


9460 

ISOURCE 


EXT 

Get_va1 ue 

! Declare 

subroutines stored 

9470 

ISOURCE 


EXT 

Get_i nf o 

! outside 

of the 

main program. 

9480 

ISOURCE 


EXT 

Get_el ement 




9490 

ISOURCE 

I 






9500 

ISOURCE 


COM 





9510 

ISOURCE 

Data_par: 

INT 

<*> 

! Declare 

c ommon 

var i ab 1 es . 

9520 

ISOURCE 

Dl_par : 

INT 

<*> 




9530 

ISOURCE 

Ns_par : 

INT 





9540 

ISOURCE 

Nt i me_par : 

INT 





9550 

ISOURCE 

! 






9560 

ISOURCE 

Arrayd: 

BSS 

39 

! Reserve 

spac e 

for data array 

9570 

ISOURCE 

El ementd: 

EQU 

Arrayd+16 

! descriptor. 

9580 

ISOURCE 

Array 1 : 

BSS 

1024 

! Reserve 

space 

for tables which hold 

9590 

ISOURCE 

Array2: 

BSS 

1024 

! the number of 

samples per count. 

9600 

ISOURCE 

Array3: 

BSS 

1024 





The subroutine works by first going to the sampl es-per-count 
tables and using them to draw the old histograms in black 
<Bit=0) to erase them. Then it calculates new samples-per- 
count tables from the data acquired from the LDV, and uses 
the new tables to draw histograms in white <Bit=l>. Then the 
subroutine returns to the main program. The very first time 
(determined by Nt i mes(=Nn> ) the subroutine is called it doesn't 
erase the old histograms because there aren't any. 


39 



9610 

ISOURCE 

Ns: 

BSS 

1 



9620 

ISOURCE 

Count : 

BSS 

1 ! Count 

and I are general purpose index uari- 

9630 

ISOURCE 

I: 

BSS 

1 ! ables 

• 

Count is usually 0, 1, or 2, to denote 

9640 

ISOURCE 

Check : 

BSS 

1 ! whether 

U,V, or W is being calculated. 

9650 

ISOURCE 

Int : 

BSS 

1 ! Int, 

Address, and Offset are all general 

9660 

ISOURCE 

Address; 

BSS 

1 ! purpose 

storage areas. 

9670 

ISOURCE 

X_coordl 

BSS 

1 



9680 

ISOURCE 

Y_coord; 

BSS 

1 



9690 

ISOURCE 

Bit: 

BSS 

1 



9700 

ISOURCE 

Nt i mes: 

BSS 

1 



9710 

ISOURCE 

j 





9720 

ISOURCE 


LIT 

30 



9730 

ISOURCE 

! 





9740 

ISOURCE 


SUB 




9750 

ISOURCE 

Draw_hi st : 

LDfi 

=Ns 

! 

Get number of samples. 

9760 

ISOURCE 


LDB 

=Ns_par 



9770 

ISOURCE 


JSM 

Get_val ue 



9780 

ISOURCE 


LDFI 

=Arrayd 

j 

Get parameters of data array. 

9790 

ISOURCE 


LDB 

= Dat a_par 



9800 

ISOURCE 


JSM 

Get_i nf o 



9810 

ISOURCE 


LDFI 

=Nt i mes 

j 

Get Ntimes, which tells if this is 

9820 

ISOURCE 


LDB 

=Nt i me_par 

I 

the first time Draw_hist is being 

9330 

ISOURCE 


JSM 

Get_oal ue 

| 

cal led. 

9840 

ISOURCE 

! 





9850 

ISOURCE 


LDFI 

= 0 

! 

If this is the first time Draw_hist 

9860 

ISOURCE 


STA 

Bi t 

j 

is being called, then jump to the 

9870 

ISOURCE 


LDA 

Nt i mes 

! 

data acquisition section. If not then 

9880 

ISOURCE 


RZA 

■*■+2 

j 

write over the old histogram with 

9890 

ISOURCE 


JMP 

Acqui re 

l 

black to erase it 

9980 

ISOURCE 

! 





9910 

ISOURCE 

Do_graph: 

LDA 

= 0 

j 

Produce histograms for all three 

9920 

ISOURCE 


STA 

Count 

| 

channel s. 

9930 

ISOURCE 

Mak e_hi st : 

LDA 

Count 

j 

Go to the end of the appropriate 

9940 

ISOURCE 


ADA 

= 1 

! 

samples-pei — count table. 

9950 

ISOURCE 


LDB 

= 1024 



9960 

ISOURCE 


STB 

I 



9970 

ISOURCE 


MPY 




9980 

ISOURCE 


ADA 

=Array 1 



9990 

ISOURCE 


ADA 

=-l 



1 0000 

ISOURCE 


STA 

Address 



10010 

ISOURCE 

Make rod: 

LDB 

Address , I 

! 

Add together each adjacent pair of 

10020 

ISOURCE 


DSZ 

Address 

! 

entries in the table. 

10030 

ISOURCE 


ADB 

Address, I 



10040 

ISOURCE 


DSZ 

Address 



10050 

ISOURCE 


DSZ 

I 



10060 

ISOURCE 


RZB 

*+2 

j 

If the result is zero, go to the next 

10070 

ISOURCE 


JMP 

Sk i p_bi ts 

! 

pair of entries. If not, make sure 

10080 

ISOURCE 


STB 

Int 

! 

the sum is <150 and then draw a 

10090 

ISOURCE 


ADB 

=-150 

j 

column with height equal to the sum. 

10100 

ISOURCE 


SBM 

»+3 



10110 

ISOURCE 


LDB 

= 150 



10120 

ISOURCE 


STB 

Int 



10130 

ISOURCE 

! 





10140 

ISOURCE 


LDA 

I 

j 

Calculate the X-coordinate of the 

10150 

ISOURCE 


SAR 

1 

! 

c o 1 umn . 

10160 

ISOURCE 


ADA 

=20 



10170 

ISOURCE 


STA 

X_coord 



10180 

ISOURCE 

! 





10190 

ISOURCE 

Mak e_bi t : 

LDA 

= 150 

j 

Calculate the Y-coordinate of the 

10200 

ISOURCE 


LDB 

Count 

1 

top of the column. 
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10210 ISOURCE 
10220 ISOURCE 
10230 ISOURCE 
10240 ISOURCE 
10250 ISOURCE 
10260 ISOURCE 
10270 ISOURCE ! 

1 0280 ISOURCE 
10290 ISOURCE 
10300 ISOURCE 
10310 ISOURCE 
10320 ISOURCE 
10330 ISOURCE 
10340 ISOURCE 
10350 ISOURCE 
10360 ISOURCE 
10370 ISOURCE 
10380 ISOURCE 
10390 ISOURCE 
10400 ISOURCE 
10410 ISOURCE 
10420 ISOURCE 
10430 ISOURCE ! 

10440 ISOURCE 
10450 ISOURCE 
10460 ISOURCE 
10470 ISOURCE 
10480 ISOURCE 
10490 ISOURCE 
10500 ISOURCE 
10510 ISOURCE 
10520 ISOURCE 
10530 ISOURCE 
10540 ISOURCE ! 

10550 ISOURCE Skip_bits: 
10560 ISOURCE 
10570 ISOURCE ! 

10580 ISOURCE 
10598 ISOURCE 
10600 ISOURCE 
10610 ISOURCE 
10620 ISOURCE 
10630 ISOURCE ! 

10640 ISOURCE Require;: 
10650 ISOURCE 
10660 ISOURCE 
10670 ISOURCE ! 

10680 ISOURCE 
10690 ISOURCE 
10700 ISOURCE Clear: 
10710 ISOURCE 
10720 ISOURCE 
10730 ISOURCE 
10740 ISOURCE ! 

10750 ISOURCE 
10760 ISOURCE 
10770 ISOURCE 
10780 ISOURCE 
10790 ISOURCE 
10800 ISOURCE 


RDB =1 
MPV 

LDB Int 
TCB 
ADR B 

STfi Y_coord 

LDfl =13 
STR Pa 
LDR =5 1 B 
SFC * 

STR R5 
LDfl Y_coord 
LDB =36 
MPY 

LDB X_coord 
SBR 4 
ADR B 
CMR 
SFC * 

STR R4 
STR R7 

LDR X_coord 
AND = 1 7B 
LDB Bit 
SBL 15 
I OR B 
SFC * 

STR R4 
STR R7 
DSZ Int 
JMP Make_bit 

DSZ I 

JMP Make_rod 

ISZ Count 
LDfl =-3 
ADR Count 
SZR *+2 
JMP Make_hist 

LDfl Bit 
SZR *+2 
JMP Stop 

LDfl =Rrray 1 
LDB =192 
CLR 16 
ADR =16 
DSZ B 
JMP Clear 

LDR =0 
STR Count 
LDfl =5 
STR Check 
LDfl Ns 
STR I 


! Prepare the graphics screen for 
! i nput . 


! Calculate word location. 


! Input word location. 


! Calculate bit location. 


! Input bit location. 


! Make histograms for all three 
! channels. 


! If Bit=l then the new histograms 
! have been drawn. Go back to the 
! main program. If not, continue. 

! Zero out the tables which hold the 
! number of samples per count. 


! Prepare to write new samples-per- 
! count tables. 
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10810 

ISOURCE 

! 



10820 

ISOURCE 

Beg i n: 

LDfl 

Count 

10830 

ISOURCE 


LDB 

Ns 

10840 

ISOURCE 


MPY 


10850 

ISOURCE 


ADA 

I 

10860 

ISOURCE 


ADA 

= -l 

10870 

ISOURCE 


STA 

E 1 ement d 

10380 

ISOURCE 


LDA 

= I nt 

10890 

ISOURCE 


LDB 

=Arrayd 

10900 

ISOURCE 


JSM 

Get_e 1 ement 

10910 

ISOURCE 


LDA 

Int 

10920 

ISOURCE 


LDB 

= 1023 

10930 

ISOURCE 


AND 

B 

10940 

ISOURCE 


LDB 

Count 

10950 

ISOURCE 


TCB 


10960 

ISOURCE 


ADB 

Check 

10970 

ISOURCE 


SZB 

Strai ght 

10980 

ISOURCE 


TCA 


10990 

ISOURCE 


ADA 

= 1024 

11000 

ISOURCE 

St rai gh " : 

STA 

Int 

11010 

ISOURCE 


LDA 

= 1024 

11020 

ISOURCE 


LDB 

Count 

11030 

ISOURCE 


MPY 


11040 

ISOURCE 


ADA 

=Array 1 

11050 

ISOURCE 


ADA 

Int 

11060 

ISOURCE 


ISZ 

A,I 

11070 

ISOURCE 

j 



11080 

ISOURCE 


ISZ 

Count 

11090 

ISOURCE 


LDA 

=-3 

11100 

ISOURCE 


ADA 

Count 

11110 

ISOURCE 


SZA 

*+2 

11120 

ISOURCE 


JMP 

Begi n 

11130 

ISOURCE 


STA 

Count 

11140 

ISOURCE 

1 



11150 

ISOURCE 


DSZ 

I 

1 1160 

ISOURCE 


JMP 

Begi n 

1 1170 

ISOURCE 

i 



11180 

ISOURCE 


ISZ 

B i t 

1 1190 

ISOURCE 


JMP 

Do graph 

1 1200 

ISOURCE 

Stop: 

RET 

1 

11210 

ISOURCE 

j 



11220 

ISOURCE 


END 

Draw hist 


! Figure out which element of the 
! data array we want to pick up. 


! Get a raw datum from the data array. 


Strip off the first six 
bits of the raw data word. 
See if Count = Check. 


If true, use the modified data 
word as an index. If not, use 
1024 minus the data word. 

Store the count we have gotten. 
Increment the appropriate place 
in the sampl es-per-count table 
by one. 


! Have U,V, and W all been done? 

! If not, go back again. 

! If so, set Count=0. 

! Continue to fill samp 1 es-per-c ount 
! tables until out of samples. 

! Now set Bit=l and go back and draw 
! the new histograms in white. 
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REM PROGRAM STflT 

! PROGRAM TO REDUCE RAW DATA FROM THREE-COMPONENT LDV SYSTEM 
! Input is 3 channels of raw LV data and calibration constants 
! Output is the three components of mean velocity, all the 
! components o*' Reynolds stress and some selected third-order 
! products. Can filter out counts which are excessively far 
! from the mean. Can display histograms of raw and 

! filtered data for all 3 channels. Plots data at end and displays 
! reduced data for entire run. Can write a summary file containing 
! all reduced data and calibration constants. 

OPTION BASE 1 
ICOM 18000 
IDELETE ALL 
I ASSEMBLE Find_vel 

COM INTEGER Data<3, 2000) , Ns, Ngs, Countbi n<3, 1024) 

COM INTEGER Range 1 , Range2 , Range3 
COM REAL Sdev,Dfl,Df2,Df3 
COM REAL Theta, Nub, Numi xl , Numi x2, Numi x3 
COM REAL Su, Sv, Sw, Suu, Svv, Sww, Suv, Suw, Svw 
COM REAL Suuv, Suw, Suuw, Suww, Suvw 
! 

DIM LI <3),R1 C3),L1 1<3>,U1 1 <3) , Seal e<3) 

DIM Date$C80 3 ,Fi 1 e$ [ 6] , Name* C 43 , T i 1 1 $ C 80] 

! 

REAL Xpos, Ypos, 2pos 
REAL Poi nt <20, 15) 

REAL Re, Ue 


INTEGER Filno,A,,Run 

j 

DATA 0,0, 0,0 

READ Max_y, Mi n_y, Fi 1 ter , Sdev 
DATA 75.64,45.64, 15.64,99. ,69. ,39. 

MAT READ LI 1 , U 1 l 
DATA 5,46,87,41,82, 123 
MAT READ LI ,R1 

DATA -.8, 2. 4,0. , . 035, -. 0025 , . 0 1 

READ LI im<l),Ul i m< 1 ) , LI i m<2) , U1 i m<2) , LI i m<3) , U1 im<3) 
DATA .4, .005,. 0025 
MAT READ Scale 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


11 ** << PROGRAM STAT : 3-COMPONENT VELOCITY DATA >> ** 
"PROGRAM STRUCTURE" 

"1. Read raw data from floppy disc" 

"2. Convert to velocity" 

"3. Calculate statistics" 

"4. Print results" 

“5. Write to disc file" 


II 


** Read raw data from floppy disc ** 


INPUT "Enter parent filename (or E to exit proqram) NameS 
IF Name$="E" THEN GOTO 1780 

INPUT "Enter no. of first and last data files ",Filel,Nf 
Ki 1 l_w=3 

INPUT "Is this a two-channel or three-channel run <2x3, Default 3) ?",Kill 


Hi st$="Y" 
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610 INPUT "Do you wish to see histograms of the data points <YxN, Default YI> ? 
" , H i st $ 

620 F i 1 ter$="N" 

630 INPUT "Do you wish to filter the output <YxN, Default N> ?", Filter* 

640 IF Filter$="N' THEN GOTO 660 

650 INPUT "Enter standard deviation to throw out " , Sdev 

660 DISP "READING RAW DATA FROM DISK" 

670 MASS STORAGE IS ":H8,0,0" 

680 PRINTER IS 0 

690 F i 1 no=F i 1 el 

700 IF F i 1 no >Nf THEN GOTO 1740 

710 Fi le$=Name$&VAL:KFi Ino) 

720 ASSIGN Fi 1 e$ TO #1 
730 ! 

740 READ #1 ; Date* 

750 READ #l;Titl$ 

760 READ #1 ; Name$ 

770 READ # 1 ; Dn 

780 READ #1 ; Nub, Num i xl , Numi x2, Numi x3, Theta, Run 

790 READ #1; Vref,Ue,Dfl,Df2,Df3 

800 READ #1 ; Xpos, Ypos, Zpos, Ns 

810 MAT READ #l;Data 

820 ! 

830 PRINT Ti 1 1 $ 

840 PRINT "DATE AND TIME OF TEST :",Date$ 

850 PRINT "FILE NAME ON DISK :",Name$ 

860 PRINT "DISK NUMBER : " , Dn 

870 PRINT "RUN NO. (POINT NO. IN PROFILE) :",Run 
880 PRINT "X,Y,Z = " , Xpos, Ypos, Zpos 
890 ! 

900 ! CALCULATE VELOCITIES 

910 DISP “CALCULATING VELOCITIES AND TAKING RUNNING SUMS" 

920 ! 

930 ICALL F i nd_ve1 
940 ! 

950 ! CALCULATE STATISTICS 

960 DISP "CALCULATING STATISTICS" 

970 ! 

980 IF Kill _w=3 THEN GOTO 1030 
990 DATA 0,0,0,0,0,13,0 
1000 RESTORE 990 

1010 READ Sw, Sww, Suw, Suw, Suuw, Suww, Suuw 
1020 ! 

1030 Ubar=SuxNgs 

1040 Vbar=SvxNgs 
1050 Wbar=Sw/Ngs 
1060 Upr i 2=SuuxNgs-Ubar*Ubar 
1070 Vpr i 2=Svv/Ngs-Vbar*Vbar 
1080 Wpr i 2=Sww/Ngs-Wbar*Wbar 
1090 Uvbar=Suv/Ngs-Ubar*Vbar 
1100 Uwbar=Suw/Ngs-Ubar*Wbar 
1110 Vwbar=SvwxNgs-Vbar*Wbar 

1120 Uuvbar=SuuV''Ngs“2*Ubar*SuvxNgs-Vbar*SuuxNgs+2*Vbar*Ubar*Ubar 

1130 Uvvbar=SuwxNgs-2* Vbar *SuvxNgs - Ubar *SvvxNgs+2*Ubar * Vbar*Vbar 

1140 Uuwbar=SuuwxNgs™2*Ubar*Suw' , Ngs-Wbar*SuuxNgs+2*Wbar*Ubar*Ubar 
1150 Uwwbar=SuwwxNgs--2*Wbar*SuwxNgs-Ubar*SwwxNgs+2*Ubar*Wbar*Wbar 

1160 Uvwbar=SuvwxNgs-Wbar*SuvxNgs-Ubar*SvwxNgs-Vbar*Suw/'Ngs+2*Ubar*Vbar*Wbar 
1170 ! 

1130 ! STORE RESULTS IN ARRAY FOR PLOTTING 

1190 ! 

1200 Fl=Fi lno-Fi lel+l 
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1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1 otr$ 
1770 
1780 
1790 
1800 


Poi nt <F1, 1 > = Ypo:» 

Poi nt <F1 , 2>=Ubar 

Poi nt <F1 , 3)=Vbar 

Poi nt <F1 , 4>=Wbar 

Point<Fl,5>=Upri2 

Point<Fl,6>=Vpri2 

Point<Fl,7)=Wpri2 

Point<Fl,8>=Uvbar 

Point<Fi,9>=Uwbar 

Poi nt <F1 , 10>=Vwbar 

Poi nt <F1 , 1 1 >=Uuvbar 

Poi nt <F1 , 12>=Uvvbar 

Poi nt <F1 , 13>=Uuwbar 

Poi nt <F1 , 14>=Uwubar 

Poi nt <F1 , 15>=Uvwbar 

IF Ypos< Max_y THEN GOTO 1390 

Max_y=Ypos 

Mi n_u=Ubar 

IF Ypos>Mi n_y THEN GOTO 1450 

M i n_y=Ypos 

Max_u=Ubar 

! PRINT RESULTS 

PRINT 

PRINT "Bragg shift frequency = " , Nub 
PRINT "Mixing f requenc i es=" , Numi xl , Numi x2, Numi x3 
PRINT "Third beam angle = ", Theta 
PRINT "Run number = " , Run 

PRINT "Tunnel reference voltage <vo1ts) = ",Vref 

PRINT "Fringe spacings = " , Df 1 , Df 2, Df 3 

PRINT "No. of samples = " , Ns 

PRINT "Ranges = ", Rangel , Range2, Range3 

PRINT "Filter Cno. of std. deviations] = ",Sdev 

PRINT “Adjusted no. of samples =",Ngs 

PRINT "Ubar = " , PROUNDCUbar , -4) 

PRINT "Vbar = " PROUND< Vbar, -4> 

PRINT "Wbar = ", PROUNDCUbar, -4) 

PRINT "Upr i 2 = 1 , PROUNDCUpr i 2, -4> 

PRINT "Vpr i 2 = " , PROUNDCVpr i 2, -4) 

PRINT " Upr i 2 = " , PROUNDCWpr i 2, -4) 

PRINT "UVbar = 1 , PROUNDCUvbar , -4) 

PRINT “UWbar = 1 , PROUNDCUwbar, -4) 

PRINT " VUbar = 1 , PROUNDCVwbar, -4) 

PRINT "UUVbar = " , PROUND < Uuvbar , -5 > 

PRINT "UWbar = " , PROUND<Uvvbar, -5> 

PRINT " UUWbar = " , PROUNDCUuwbar, -5> 

PRINT "UWWbar = " , PROUNDCUwwbar, -5) 

PRINT "UVUbar = " , PROUNDCUvwbar, -5> 

PRINT 

IF Hist$="Y" THEN GOSUB Histogram 

F i 1 no=F i 1 no+1 

GOTO 700 

Ue=Max_u-M i n_u 

Plotr$="Y" 

INPUT " Do you wish to see the turbulence quantities <YxN, Default Y> ?",P 

IF P1otr$="Y" THEN GOSUB Plot 
MASS STORAGE IS ":H8,0,1" 

PRINTER IS 16 
END 
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1310 ! 

1820 ! **********************PLOT DATA***************************************** 

1830 ! 

1840 Plot: PRINTER IS 16 

1350 PRINT "Estimate of normalizing velocity is : " ; Ue 

I860 INPUT "Enter correct normalizing velocity if d i f f erent . " , Ue 

1870 Graph$="Y" 

1880 INPUT "Do you wish to graph the turbulence quantities (Y/N, Default Y) ?", 
Graph* 

1890 PRINTER IS 0 

1900 PRINT "Normalizing velocity is : " ; Ue 
1910 ! 

1920 FOR 1=1 TO 4 
1930 DATA 1,2, 2, 3 
1940 READ Exp 
1950 Kn=3 

I960 IF 1=4 THEN Kn=5 
1970 FOR J=1 TO Nf-Nfile+1 
1980 FOR K= 1 TO Kn 
1990 Index=<I-l)*3+K+l 

2008 Poi nt ( J, Index)=Poi nt < J, IndexVUe-'Exp 

2010 IF Index>10 THEN Poi nt ( J, Index)=Poi nt < J, Index>*1000 

2020 NEXT K 

2030 NEXT J 

2040 NEXT I 

2050 IF Graph*="N" THEN GOTO 2860 
2060 ! 

2070 PLOTTER IS 13, "GRAPHICS" 

2080 GRAPHICS 
2090 DEG 
2100 ! 

2110 Ymax=Poi nt < 1 , 1 ) 

2120 Ymin=Ymax 

2130 FOR J=1 TO Nf-Filel+1 

2140 IF Poi nt < J, 1 ) >Ymax THEN Ymax=Po i nt < J , 1 > 

2150 IF Poi nt < J, 1 >< Ymi n THEN Ym i n=Po i nt < J , 1 ) 

2160 NEXT J 
2170 ! 

2130 FOR 1=1 TO 3 
2190 LIMIT 0,184.47,0,149.8 
2200 LOCATE LI <I),R1 < I >,39, 99 
2210 CLIP LI < I ) , R1 < I > , 39, 99 
2220 ! 

2230 Xmax=Point<l, <I-l)*3+2> 

2240 Xmin=Xmax 

2250 FOR J=1 TO Nf-Filel+1 

2260 FOR K= 1 TO 3 

2270 IF Poi nt < J, < 1-1 >*3+K+l ) >Xmax THEN Xmax=Poi nt ( J, C I - 1 >*3+K+l > 

2280 IF Poi nt < J, < 1-1 )*3+K+l)<Xmi n THEN Xrai n=Poi nt < J, ( 1-1 >*3+K+l > 

2290 NEXT K 
2300 NEXT J 
2310 ! 

2320 SCALE L 1 i m < I > , U 1 i m < I) , Ym i n , Ymax 
2330 FRAME 

2340 AXES Scaled), . 1,0,0 
2350 UNCLIP 
2360 ! 

2370 LORG 3 
2380 CSIZE 2.5 
2390 Ypos=Ymin 

2400 MOVE LI im<I)+Scale(I)x5, Ypos 
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2410 DRAW LI i m< I ) , Ypos 
2420 SETGU 
2430 RPLOT 2, 0,-2 
2440 SETUU 

2450 LABEL DROUND< Ypos, 2) 

2460 Ypos=Ypos+. 1 

2470 IF Ypos<=Yma.x THEN GOTO 2400 
2480 ! 

2490 LORG 6 
2500 Xpos=Ll i m< I > 

2510 MOVE Xpos, Ymin-,04 
2520 LABEL DROUNDCXpos, 2) 

2530 MOVE Xpos , Ymi n 
2540 DRAW Xpos , Ym i n+ . 05 
2550 Xpos=Xpos+Scal e< I ) 

2560 IF Xpos<Ul i m< I ) THEN GOTO 2510 
2570 ! 

2580 FOR J=1 TO Nf-Filel+1 
2590 Ypos = Poi nt < J, 1 ) 

2600 FOR K= 1 TO 3 

2610 MOVE PointCJ, <I-1)*3+K+1), Ypos 
2620 SETGU 

2630 ON K GOTO 2640,2700,2760 

2640 FOR Arc =0 TO 360 STEP 20 

2650 PDIR Arc 

2660 RPLOT .5,0 

2670 NEXT Arc 

2680 PDIR 0 

2690 GOTO 2810 

2700 RPLOT .5,. 5, -2 

2710 RPLOT .5, -.5,-1 

2720 RPLOT -.5, -.5,-1 

2730 RPLOT -.5, .5,-1 

2740 RPLOT .5,. 5,-1 

2750 GOTO 2810 

2760 RPLOT 0, .5,-2 

2770 RPLOT 0, -. 5, -1 

2780 RPLOT .5, 0,-2 

2790 RRLOT -.5, 0,-1 

2800 GOTO 2810 

2810 SETUU 

2820 NEXT K 

2830 NEXT J 

2840 NEXT I 

2850 DUMP GRAPHICS 

2860 IMAGE 2X"Y"10X, 'u2/Uo2",3X, "o2/Uo2" , 3X, "w2xUo2",5X, "uv/'Uo2“ , 3X, "uw/Uo2",3X 
, ,, ou)/Uo2",/ 

2870 PRINT USING 2860 
2880 FOR 1=1 TO Nf-Filel+1 

2890 IMAGE MDD. DD, 6X, MD. DDDD, 2X, MD. DDDD, 2X, MD. DDDD, 4X, MD. DDDD, 2X, MD. DDDD, 2X, MD. 
DDDD 

2900 PRINT USING 2890; Poi nt < I , 1 ) , Po i nt < I , 5) , Po i nt < I , 6 ) , Poi nt ( I , 7) , Po i nt < I , 8) , Po 
int<I,9),Point<I, 10) 

2910 NEXT I 

2920 IMAGE 2/,3X,"Y", 9X, "uuv/Uo3“ , 3X, "uw/Uo3" , 5X, "uuwxUo3" , 3X, "uuiw/'Uo3" , 5X, "u 
v w x U o 3 " , x 

2930 PRINT USING 2920 
2940 FOR 1=1 TO Nf-File+1 

2950 IMAGE MDD. DD, 6X : , MD. 5D, 2X, MD. 5D, 4X, MD. 5D, 2X, MD. 5D, 4X, MD. 5D 

2960 PRINT USING 2950; Poi nt < 1 , 1 ) , Poi nt < 1 , 1 1 ) , Poi nt < 1 , 12) , Poi nt <1 , 13) , Poi nt < 1 , 14 
) , Po i nt < 1 , 15) 
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NEXT I 


2970 
2980 

2990 ! WRITE SUMMARY DATA FILE 

3000 ! 

3010 Sum$="N" 

3020 INPUT "Do you wish to write a Summary Data File <YxN Default N> ?",Sum$ 
3030 IF Sum$="N" THEN GOTO 3210 . 

3040 Fi1e$=Name$ 

3050 DISP "File ";File$;" being written to disk" 

3060 MASS STORAGE IS ":H8,0,0" 

3070 CREATE Fi 1 e*, 40 
3080 ASSIGN Fi 1e$ TO #1 
3090 PRINT #1 ; DateJ 
3100 PRINT #l;Titl$ 

3110 PRINT #1 ; Name* 

3120 PRINT # 1 ; Dn 

3130 PRINT #1 ; Nub, Numi xl , Numi x2, Numi x3, Theta, Run 

3140 PRINT #l;Vref,Ue,Dfl,Df2,Df3 

3150 PRINT #1 ; Ns, Fi 1 el , Nf 

3160 FOR 1=1 TO Nf-Nfile+1 

3170 FOR J=1 TO 15 

3180 PRINT #1 ; Poi nt < I , J) 

3190 NEXT J 
3200 NEXT I 
3210 RETURN 
3220 ! 

3230 ! *********************DRAW HISTOGRAMS************************************ 

3240 ! 

3250 Histogram: PLOTTER IS 13, "GRAPHICS" 

3260 GRAPHICS 
3270 K= 1 

3280 FOR 1=1 TO 3 
3290 Cmax=0 

3300 FOR J=1 TO 1023 STEP 2 

3310 Count bi n< I , J)=Countbi n < I , J)+Countbin< I , J+l ) 

3320 IF Countbi n< I , J > >Cmax THEN Cmax=Countbi n< I , J) 

3330 NEXT J 

3340 IF <Cmax=0> OR <Cmax=2000> THEN GOTO 3660 
3350 LIMIT 0,184.47,0,149.8 
3360 LOCATE 7. 05, 1 19 . 63, LI 1 CIO , U1 1 <IO 
3370 CLIP 7.05, 119. 63, LI 100,111 1 < K> 

3380 K=K+ 1 

3390 SCALE 1 , 1 024 , 0 , Cmax 

3400 FRAME 

3410 UNCLIP 

3420 LINE TYPE 1 

3430 CSI2E 2.8 

3440 LORG 3 

3450 FOR Y= 1 TO 4 

3460 Ypos=Y*Cmaxx4 

3470 PLOT 1 0 , Ypos , -2 

3480 PLOT -9, Ypos, -1 

3490 MOVE 3, Ypos 

3500 LABEL PROUND < Ypos , 0 ) 

3510 NEXT Y 

3520 LORG 6 

3530 Ypos=Cmax/'20 

3540 FOR X=1 TO 8 

3550 Xpos=X*128 

3560 PLOT Xpos, Ypos , -2 
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3570 PLOT Xpos , -Ypos , - 1 
3580 MOVE Xpos , -Ypos*l . 3 
3590 LABEL Xpos 
3600 NEXT X 

3610 FOR J=1 TO 1023 STEP 2 

3620 IF CountbinCI, J>=0 THEN GOTO 3650 

3630 MOVE J,0 

3640 DRAW J, Countbi nil , J) 

3650 NEXT J 
3660 NEXT I 
3670 MOVE 512, -Ypos*5 
3680 LORG 6 

3690 LABEL "PLOT OF COUNT VS NUMBER OF SAMPLES PER COUNT" 

3700 DUMP GRAPHICS -Ypos*10 
3710 RETURN 
3720 ! 

3730 ! 

3740 ISOURCE NAM Find_vel 

3750 ISOURCE ! 

3760 ISOURCE ! This subroutine converts raw data counts to instantaneous 
3770 ISOURCE ! velocities, then sums several different products of the 
3780 ISOURCE ! velocity components. All input and output data is passed 
3790 ISOURCE ! through the COMMON stoarage area. The inputs are the raw 
3800 ISOURCE ! data array (Arrayd), the Bragg shift frequency <Nub), the 
3810 ISOURCE ! mixing frequencies (Numixl, Numix2, Numix3>, the fringe 
3820 ISOURCE ! spacings <Dfl, Df2, Df3>, the crossing angle of the third 
3830 ISOURCE ! beam ITheta), and the number of samples in a data point <Ns>. 
3840 ISOURCE ! The outputs are the summations of various products of the 
3850 ISOURCE ! velocity components, including U, V, W, U*U, V*V, W*W, 

3860 ISOURCE ! U*V, IJ*W, V*W, U*U*V, U*V*V, U*U*W, U*W*W, and U*V*W. 

3870 ISOURCE ! 

3880 ISOURCE EXT Get_value ! Declare subroutines stored 

3890 ISOURCE EXT Get_i nf o ! outside of the main program. 

3900 ISOURCE EXT Get_e1ement 

3910 ISOURCE EXT Put_element 

3920 ISOURCE EXT Int_to_rel 

3930 ISOURCE EXT Rel_to_int 

3940 ISOURCE EXT Rel_math 

3950 ISOURCE EXT Put_value 

3960 ISOURCE ! 

3970 ISOURCE COM 


3980 

ISOURCE 

Dat a_par : 

INT 

<*> 

3990 

ISOURCE 

Ns_par : 

INT 


4000 

ISOURCE 

Ngs_par :: 

INT 


4010 

ISOURCE 

Cbi n_par : 

INT 

<*> 

4020 

ISOURCE 

Ranl_par: 

INT 


4030 

ISOURCE 

Ran2_p ar 1 

INT 


4040 

ISOURCE 

Ran3_par : 

INT 


4050 

ISOURCE 

Sdev_par : 

REL 


4060 

ISOURCE 

Df l_par 5 

REL 


4070 

ISOURCE 

Df2_par 3 

REL 


4080 

ISOURCE 

Df3_par ;! 

REL 


4090 

ISOURCE 

Thet a_par : 

REL 


4100 

ISOURCE 

Nub_par ; 

REL 


4110 

ISOURCE 

Nmi xl_par : 

REL 


4120 

ISOURCE 

Nmi x2_par : 

REL 


4130 

ISOURCE 

Nmi x3_par : 

REL 


4140 

ISOURCE 

Su_par : 

REL 


4150 

ISOURCE 

Sv_par : 

REL 


4160 

ISOURCE 

Sw_par : 

REL 



Declare common variables. 
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4170 

ISOURCE 

Suu_par : 

REL 



4180 

ISOURCE 

Swv_parl 

REL 



4190 

ISOURCE 

Suu_par 

REL 



4200 

ISOURCE 

Suu_par 1 

REL 



4210 

ISOURCE 

Suu_par 5 

REL 



4220 

ISOURCE 

Suu_par : 

REL 



4230 

ISOURCE 

Suuu_par : 

REL 



4240 

ISOURCE 

Suuu_par : 

REL 



4250 

ISOURCE 

Suuu par: 

REL 



4260 

ISOURCE 

Suuu_par : 

REL 



4270 

ISOURCE 

Suuu_par : 

REL 



4280 

ISOURCE 

! 




4290 

ISOURCE 

Arrayd: 

BSS 

39 

! Reserue space for data array 

4300 

ISOURCE 

El ementd: 

EQU 

Arrayd+16 ! descriptor. 

4310 

ISOURCE 

Array 1 : 

BSS 

4096 

! Reserue space for lookup tables 

4320 

ISOURCE 

Array2 : 

BSS 

4096 

! for count to uelocity couersion. 

4330 

ISOURCE 

Array3: 

BSS 

4096 


4340 

ISOURCE 

B i n_u : 

BSS 

1024 

! Reserue space for table of counts 

4350 

ISOURCE 

B i n_u : 

BSS 

1024 

! used to generate histograms. 

4360 

ISOURCE 

B i n_u : 

BSS 

1024 


4370 

ISOURCE 

Rangel : 

BSS 

1 

! Reserue space for uarious input 

4380 

ISOURCE 

Range2: 

BSS 

1 

! and output uariables. 

4390 

ISOURCE 

Range3: 

BSS 

1 


4400 

ISOURCE 

Ns: 

BSS 

1 


4410 

ISOURCE 

Df l: 

BSS 

4 


4420 

ISOURCE 

Df2: 

BSS 

4 


4430 

ISOURCE 

Df3: 

BSS 

4 


4440 

ISOURCE 

Thet a: 

BSS 

4 


4450 

ISOURCE 

Nub: 

BSS 

4 


4460 

ISOURCE 

Numi xl : 

BSS 

4 


4470 

ISOURCE 

Nunn x2: 

BSS 

4 


4480 

ISOURCE 

Numi x3: 

BSS 

4 


4490 

ISOURCE 

Su: 

BSS 

4 


4500 

ISOURCE 

Su: 

BSS 

4 


4510 

ISOURCE 

Sw: 

BSS 

4 


4520 

ISOURCE 

Suu : 

BSS 

4 


4530 

ISOURCE 

Suu: 

BSS 

4 


4540 

ISOURCE 

Suu: 

BSS 

4 


4550 

ISOURCE 

Suu: 

BSS 

4 


4560 

ISOURCE 

Suu: 

BSS 

4 


4570 

ISOURCE 

Suu: 

BSS 

4 


4580 

ISOURCE 

Suuu: 

BSS 

4 


4590 

ISOURCE 

Suuu : 

BSS 

4 


4600 

ISOURCE 

Suuu : 

BSS 

4 


4610 

ISOURCE 

Suuu: 

BSS 

4 


4620 

ISOURCE 

Suuu: 

BSS 

4 


4630 

ISOURCE 

Count : 

BSS 

1 ! 

Count and I are general purpose index uari- 

4640 

ISOURCE 

i: 

BSS 

1 ! 

ables. Count is usually 0, 1, or 2 to denote 

4650 

ISOURCE 

Check : 

BSS 

1 ! 

uhether U, V, or W is being calculated. 

4660 

ISOURCE 

Int : 

BSS 

1 ! 

Int, Address, and Offset are all general 

4670 

ISOURCE 

I nt 2 : 

BSS 

1 ! 

purpose storage areas. 

4680 

ISOURCE 

Address: 

BSS 

1 


4690 

ISOURCE 

Of fset : 

BSS 

1 


4700 

ISOURCE 

Rl: 

BSS 

4 ! 

Rl, R2, and R3 are the count-to- 

4710 

ISOURCE 

R2: 

BSS 

4 ! 

frequency conuersion factors. 

4720 

ISOURCE 

R3: 

BSS 

4 


4730 

ISOURCE 

Xuar : 

BSS 

4 ! 

Xuar and Yuar are general purpose real 

4740 

ISOURCE 

Yuar : 

BSS 

4 ! 

number storage areas. 

4750 

ISOURCE 

u: 

BSS 

4 ! 

U, V, and W are the instantaneous uelocity 

4760 

ISOURCE 

v: 

BSS 

4 ! 

c omponent s . 
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4770 

ISOURCE 

w: 

BSS 

4 

4780 

ISOURCE 

Uu: 

BSS 

4 

4790 

ISOURCE 

Uv: 

BSS 

4 

4800 

ISOURCE 

Cos: 

BSS 

4 ! Cos and Sin are the cos and sin of Theta 

4810 

ISOURCE 

Sin: 

BSS 

4 

4820 

ISOURCE 

Max_u: 

BSS 

1 

4830 

ISOURCE 

Max_v: 

BSS 

1 

4340 

ISOURCE 

Max_w: 

BSS 

1 

4850 

ISOURCE 

Rvg_u: 

BSS 

4 

4860 

ISOURCE 

flvg_v: 

BSS 

4 

4870 

ISOURCE 

flvg_u: 

BSS 

4 

4880 

ISOURCE 

Dev_u! 

BSS 

4 

4890 

ISOURCE 

Dev_v! 

BSS 

4 

4900 

ISOURCE 

Dev_w: 

BSS 

4 

4910 

ISOURCE 

Ngs: 

BSS 

1 

4920 

ISOURCE 

Sdev : 

BSS 

4 

4930 

ISOURCE 

Rad: 

CRT 

5. 729578E 1 

4940 

ISOURCE 

Mill: 

DAT 

1.E6 

4950 

ISOURCE 

One: 

DAT 

1 . 

4960 

ISOURCE 

Zero: 

DAT 

0. 

4970 

ISOURCE 


LIT 

80 

4980 

ISOURCE 

! 



4990 

ISOURCE 


SUB 


5000 

ISOURCE 

F i nd_ve 1 : 

LDA 

=Ns ! Get number of samples. 

5010 

ISOURCE 


LDB 

=Ns_par 

5020 

ISOURCE 


JSM 

Get_val ue 

5030 

ISOURCE 


LDR 

=flrrayd ! Get parameters of data array. 

5040 

ISOURCE 


LDB 

=Data_par 

5050 

ISOURCE 


JSM 

Get_info 

5060 

ISOURCE 


LDA 

=Sdev ! Get input parameters. 

5070 

ISOURCE 


LDB 

=Sdev_par 

5080 

ISOURCE 


JSM 

Get value 

5090 

ISOURCE 


LDA 

=Df 1 

5100 

ISOURCE 


LDB 

=Dfl_par 

51 10 

ISOURCE 


JSM 

Get value 

5120 

ISOURCE 


LDR 

=Df 2 

5130 

ISOURCE 


LDB 

=Df 2_par 

5140 

ISOURCE 


JSM 

Get value 

5150 

ISOURCE 


LDR 

= Df 3 

5160 

ISOURCE 


LDB 

=Df 3_par 

5170 

ISOURCE 


JSM 

Get_value 

5180 

ISOURCE 


LDR 

=Thet a 

5190 

ISOURCE 


LDB 

=Thet a_par 

5200 

ISOURCE 


JSM 

Get_va1 ue 

5210 

ISOURCE 


LDR 

= Nub 

5220 

ISOURCE 


LDB 

=Nub_par 

5230 

ISOURCE 


JSM 

Get_val ue 

5240 

ISOURCE 


LDR 

=Numi xl 

5250 

ISOURCE 


LDB 

=Nm i xi_par 

5260 

ISOURCE 


JSM 

Get_val ue 

5270 

ISOURCE 


LDR 

=Num i x2 

5280 

ISOURCE 


LDB 

=Nmi x2_par 

5290 

ISOURCE 


JSM 

Get_value 

5300 

ISOURCE 


LDR 

=Numi x3 

5310 

ISOURCE 


LDB 

=Nmi x3_par 

5320 

ISOURCE 


JSM 

Get_val ue 

5330 

ISOURCE 

j 



5340 

ISOURCE 

! The 1 oop 

headed by Get freq is repeated three times to get 

5350 

ISOURCE 

! the count 

-to- 

■frequency conversion factors <which depend on 

5360 

ISOURCE 

! the range) for U,V, and W. Whenever a loop is controlled by 


51 



5370 

ISOURCE 

! the variable 

"Count", th 

5380 

5390 

5400 

ISOURCE 

ISOURCE 

ISOURCE 

! the same 
! 

! 

for 

U,V, and W. 

5410 

ISOURCE 


LDR 

= 0 

5420 

ISOURCE 


STfl 

Count 

5430 

ISOURCE 

Get_f req : 

LDfl 

Count 

5440 

5450 

ISOURCE 

ISOURCE 


LDB 

MPY 

Ns 

5460 

ISOURCE 


STfl 

El ementd 

5470 

ISOURCE 


LDfl 

= I nt 

5480 

ISOURCE 


LDB 

=flrrayd 

5490 

ISOURCE 


JSM 

Get_e 1 ement 

5508 

ISOURCE 


LDfl 

I nt 

5510 

ISOURCE 


LDB 

=15360 

5520 

ISOURCE 


AND 

B 

5530 

5540 

ISOURCE 

ISOURCE 


SflR 

TCR 

10 

5550 

ISOURCE 


LDB 

= 15 

5560 

ISOURCE 


ADA 

B 

5570 

ISOURCE 


LDB 

=Rangel 

5580 

ISOURCE 


ADB 

Count 

5590 

ISOURCE 


STfl 

B, I 

5600 

ISOURCE 


LDB 

= 1 

5610 

ISOURCE 


SZfl 

Loopend 

5620 

ISOURCE 

Loop: 

SBL 

1 

5630 

ISOURCE 


DSZ 

fi 

5640 

ISOURCE 


JMP 

Loop 

5650 

ISOURCE 

Loopendl 

STB 

Int 

5660 

ISOURCE 


LDfl 

= I nt 

5670 

ISOURCE 


STfl 

Oper_l 

5680 

ISOURCE 


LDfl 

=Yvar 

5690 

ISOURCE 


STfl 

Resul t 

5700 

ISOURCE 


JSM 

Int to rel 

5710 

ISOURCE 


LDfl 

==3. 2E4 

5720 

ISOURCE 


LDB 

=Xvar 

5730 

ISOURCE 


XFR 

4 

5740 

ISOURCE 


STB 

Oper_l 

5750 

ISOURCE 


LDfl 

=Yvar 

5760 

ISOURCE 


STfl 

0per_2 

5770 

ISOURCE 


LDfl 

Count 

5780 

ISOURCE 


SflL 

2 

5790 

ISOURCE 


ADA 

= R 1 

5800 

ISOURCE 


STfl 

Resul t 

5810 

ISOURCE 


LDfl 

=2 

5820 

ISOURCE 


LDB 

=147155B 

5830 

ISOURCE 


JSM 

Rel_math 

5840 

ISOURCE 


ISZ 

Count 

5850 

ISOURCE 


LDfl 

=3 

5860 

ISOURCE 


CPfl 

Count 

5870 

ISOURCE 


JMP 

*+2 

5880 

5890 

ISOURCE 

ISOURCE 

! 

JMP 

Get_f req 

5900 

ISOURCE 


LDfl 

=flrray 1 

5910 

ISOURCE 


LDB 

=768 

5920 

ISOURCE 

Cont i nue: 

CLR 

16 

5930 

ISOURCE 


ADA 

= 16 

5940 

ISOURCE 


DSZ 

B 

5950 

5960 

ISOURCE 

ISOURCE 

, 

JMP 

Cont i nue 


e loop contains operations which are 


Get the first word of the column of 
the data array which contains the 
velocity component for which we want 
to get the range. 


! Mask and rotate to get the four 
! bits containing the range. 


! Subtract from 15 to get the 
! actual range. 

! Store the actual range in RangeN 
! so that we can transfer it to the 
! main program 

! Use the range to find the power 
! of two needed for the divisor. 


! Convert the power of two into a 
! real number. 


! Divide 3.2E4 by the appropriate 
! power of two, using BCD math. 


! Decide whether to put the result 
! in Rl, R2, or R3, depending on Count 


! Now, finally, call the utility to 
! perform the division. 

! Increment and check Count so as 
! to follow the loop three times. 


Zero out the entire count-to- 
velocity conversion table so that 
it must be recalculated for each 
point. (This must be done if the 
mixing frequency or ranges are 
changed between counts.) 
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5970 

ISOURCE 


LDR 

= Su 

! Set initial values of Su, Sv, 

5930 

ISOURCE 


LDB 

= 15 

! Sw, Suu, etc. to zero. 

5990 

ISOURCE 

Clear: CLR 

4 


6G00 

ISOURCE 


ADR 

=4 


6010 

ISOURCE 


DS2 

B 


6020 

ISOURCE 


JMP 

Clear 


6030 

ISOURCE 

j 




6040 

ISOURCE 


LDR 

=Bi n u 

! Clear the areas of memory which 

6050 

ISOURCE 


LDB 

= 192 

! will be used to hold the number of 

6060 

ISOURCE 

Cl ear2: CLR 

16 

! samples per count for the histogram 

6070 

ISOURCE 


RDfl 

= 16 


6080 

ISOURCE 


DSZ 

B 


6090 

ISOURCE 


JMP 

Clear2 


6100 

ISOURCE 

j 




6110 

ISOURCE 


LDR 

= Thet a 

! Convert Theta from degrees to 

6120 

ISOURCE 


STR 

Oper 1 

! radians using the Rel_math 

6130 

ISOURCE 


LDR 

=Rad 

! utility. 

6140 

ISOURCE 


STR 

0per_2 


6150 

ISOURCE 


LDR 

=Xvar 


6160 

ISOURCE 


STR 

Resul t 


6170 

ISOURCE 


LDR 

=2 


6180 

ISOURCE 


LDB 

=1471 55B 


6190 

ISOURCE 


JSM 

Rel_math 


6200 

ISOURCE 

j 




6210 

ISOURCE 


LDR 

=Xvar 

! Find the sine and cosine of 

6220 

ISOURCE 


STR 

Oper_l 

! Theta, and store them in. the 

6230 

ISOURCE 


LDR 

=Si n 

! locations Sin and Cos, respectively 

6240 

ISOURCE 


STR 

Resul t 


6250 

ISOURCE 


LDR 

= 1 


6260 

ISOURCE 


LDB 

=34213B 


6270 

ISOURCE 


JSM 

Rel_math 


6280 

ISOURCE 


LDR 

=Cos 


6290 

ISOURCE 


STR 

Resu 1 1 


6300 

ISOURCE 


LDR 

= 1 


6310 

ISOURCE 


LDB 

=34224B 


6320 

ISOURCE 


JSM 

Rel_math 


6330 

ISOURCE 

l 




6340 

ISOURCE 

j 

The loop defi 

ned by Get 

int calculates intermediate values 

6350 

ISOURCE 

! 

used in converting counts to velocities. These values are the 

6360 

ISOURCE 

| 

same -'or al 1 

samples in 

a point, so they can be calculated 

6370 

ISOURCE 

j 

separately. The loop calculates (Nub-NumixN) and <Mill*DfN>, 

6380 

ISOURCE 

j 

where N is 1, 

2, and 3. 


6390 

ISOURCE 


LDR 

=3 


6400 

ISOURCE 


STR 

Count 


6410 

ISOURCE 

Get int;i LDR 

Count 


6420 

ISOURCE 


RDR 

=-l 


6430 

ISOURCE 


SRL 

2 


6440 

ISOURCE 


STR 

Offset 


6450 

ISOURCE 


LDR 

=Nub 

! Find Numi xN=Nub-Numi xN. 

6460 

ISOURCE 


STR 

Oper_l 


6470 

ISOURCE 


LDR' 

=Numi xl 


6480 

ISOURCE 


RDR 

Offset 


6490 

ISOURCE 


STR 

0per_2 


6500 

ISOURCE 


STR 

Resul t 


6510 

ISOURCE 


LDR 

=2 


6520 

ISOURCE 


LDB 

=146717B 


6530 

ISOURCE 


JSM 

Rel math 


6540 

ISOURCE 


LDR 

=Df 1 

! Find Df N=Mi 1 l*DfN. 

6550 

ISOURCE 


RDR 

Offset 


6560 

ISOURCE 


STR 

Oper_l 
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6570 

ISOURCE 


STfl 

Resul t 



6580 

ISOURCE 


LDfl 

=M i 1 1 



6590 

ISOURCE 


STfl 

Oper 2 



6600 

ISOURCE 


LDfl 

=2 



6610 

ISOURCE 


LDB 

= 1 47037B 



6620 

ISOURCE 


JSM 

Re1_math 



6630 

ISOURCE 


DSZ 

Count 



6640 

ISOURCE 


JMP 

Get_i nt 



6650 

ISOURCE 

J 





6660 

ISOURCE 

j 

This section 

f i 1 ters the 

data to remove counts which are 

6670 

ISOURCE 

! 

excessi vel y 

far from the 

mean. First the section creates 

6600 

ISOURCE 

! 

a list of the number of 

times each particular count appears 

6690 

ISOURCE 

! 

The list consists of 3 1024 word arrays (called Bin_u, Bin_v, 

6700 

ISOURCE 

| 

and Bin w ) where the address of each word corresponds to the 

6710 

ISOURCE 

j 

count, and the value of 

the word indicates the number of times 

6720 

ISOURCE 

j 

that particular count has 

appeared. 

6730 

ISOURCE 

i 





6740 

ISOURCE 


LDfl 

= 0 



6750 

ISOURCE 


STfl 

Count 



6760 

ISOURCE 


LDfl 

= 5 

J 

Set Check to determine whether the 

6770 

ISOURCE 


STfl 

Check 

1 

count or 1024 minus the count is to 

6780 

ISOURCE 


LDfl 

Ns 

I 

be used. 

6790 

ISOURCE 


STfl 

I 



6300 

ISOURCE 

! 





6810 

ISOURCE 

Fi 

11 bin: LDfl 

Count 



6820 

ISOURCE 


LDB 

Ns 

! 

Figure out which element of the 

6830 

ISOURCE 


MPY 


j 

data array we want to pick up. 

6840 

ISOURCE 


ADA 

I 



6850 

ISOURCE 


ADA 

=-l 



6860 

ISOURCE 


STfl 

E 1 ement d 



6870 

ISOURCE 


LDfl 

= I nt 

! 

Get a raw datum from the data array. 

6880 

ISOURCE 


LDB 

=flrrayd 



6890 

ISOURCE 


JSM 

Get_el ement 



6900 

ISOURCE 


LDfl 

Int 



6910 

ISOURCE 


LDB 

= 1023 

j 

Strip off the first six 

6920 

ISOURCE 


AND 

B 

i 

bits of the raw data word. 

6930 

ISOURCE 


LDB 

Count 

j 

See if Count = Check. 

6940 

ISOURCE 


TCB 




6950 

ISOURCE 


flDB 

Check 



6960 

ISOURCE 


S2B 

Sk i p 

! 

If true, use the modified data 

6970 

ISOURCE 


TCfl 


j 

word as an index. If not, use 

6980 

ISOURCE 


ADA 

= 1024 

! 

1024 minus the data word. 

6990 

ISOURCE 

Skip: STfl 

Int 

j 

Store the count we have gotten. 

7000 

ISOURCE 


LDfl 

Count 



7010 

ISOURCE 


LDB 

= 1024 



7020 

ISOURCE 


MPY 




7030 

ISOURCE 


ADA 

Int 

! 

Go to the address in the array 

7040 

ISOURCE 


ADA 

= -l 

j 

correspondi ng to the value of 

7050 

ISOURCE 


ADA 

=Bi n u 

j 

the count. 

7060 

ISOURCE 


IS2 

A, I ~ 

! 

Increment the word by one, 

7070 

ISOURCE 


IS2 

Count 

j 

indicating that one more count 

7080 

ISOURCE 


LDfl 

=-3 

1 

with that value has been read. 

7090 

ISOURCE 


ADA 

Count 



7100 

ISOURCE 


SZfl 

*+2 



7110 

ISOURCE 


JMP 

Fill b i n 



7120 

ISOURCE 


STfl 

Count 



7130 

ISOURCE 


DS2 

I 



7140 

ISOURCE 


JMP 

F i 1 1 _b i n 



7150 

ISOURCE 

! 





7160 

ISOURCE 

! 

When Sdev>0, 

use this section to filter the data. 
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7170 

ISOURCE 

! 

The mean 

and 

standard de 

7180 

ISOURCE 

j 

The program 

then throws 

7190 

ISOURCE 

| 

Sdev standard deviations 

7200 

ISOURCE 

! 




7210 

ISOURCE 

Check dev: 

LDB 

= S d e v 

7220 

ISOURCE 



RDB 

= 1 

7230 

ISOURCE 



LDR 

B, I 

7240 

ISOURCE 



RZfl 

*+2 

7250 

ISOURCE 



JMP 

Set_begi n 

7260 

ISOURCE 

j 




7270 

ISOURCE 



LDR 

= 0 

7280 

ISOURCE 



STR 

Count 

7290 

ISOURCE 

Get avg;! 

LDR 

=Zero 

7300 

ISOURCE 



LDB 

Count 

7310 

ISOURCE 



SBL 

2 

7320 

ISOURCE 



RDB 

=Rvg u 

7330 

ISOURCE 



XFR 

4 

7340 

ISOURCE 



LDR 

Count 

7350 

ISOURCE 



ADR 

= 1 

7360 

ISOURCE 



LDB 

= 1024 

7370 

ISOURCE 



MPY 


7380 

ISOURCE 



RDR 

=Bi n_u 

7390 

ISOURCE 



STR 

Address 

7400 

ISOURCE 



LDR 

= 1023 

7410 

ISOURCE 



STR 

I 

7420 

ISOURCE 

Sum count: 

LDR 

= 1 

7430 

ISOURCE 



STR 

Oper_l 

7440 

ISOURCE 



LDR 

=Xvar 

7450 

ISOURCE 



STR 

Resul t 

7460 

ISOURCE 



JSM 

Int_t o_rel 

7470 

ISOURCE 

j 




7480 

ISOURCE 



LDR 

Address 

7490 

ISOURCE 



STR 

Oper_l 

7500 

ISOURCE 



LDR 

=Yvar 

7510 

ISOURCE 



STR 

Resu 1 t 

7520 

ISOURCE 



JSM 

Int_t o_rel 

7530 

ISOURCE 

j 




7540 

ISOURCE 



LDR 

=Xvar 

7550 

ISOURCE 



STR 

Oper_l 

7560 

ISOURCE 



LDR 

=Yvar 

7570 

ISOURCE 



STR 

0per_2 

7580 

ISOURCE 



STR 

Resul t 

7590 

ISOURCE 



LDR 

=2 

7600 

ISOURCE 



LDB 

=147037B 

7610 

ISOURCE 



JSM 

Rel_math 

7620 

ISOURCE 

j 




7630 

ISOURCE 



LDB 

Count 

7640 

ISOURCE 



SBL 

2 

7650 

ISOURCE 



RDB 

=Rvg_u 

7660 

ISOURCE 



STB 

Oper_l 

7670 

ISOURCE 



STB 

Resu 1 t 

7680 

ISOURCE 



LDR 

=Yvar 

7690 

ISOURCE 



STR 

Oper 2 

7700 

ISOURCE 



LDR 

=2 

7710 

ISOURCE 



LDB 

=146721B 

7720 

ISOURCE 



JSM 

Rel_math 

7730 

ISOURCE 

! 




7740 

ISOURCE 



DSZ 

Address 

7750 

ISOURCE 



DSZ 

I 

7760 

ISOURCE 



JMP 

Sum count 


viation of the data are calculated, 
out all the data which are farther than 
from the mean. 


If Sdev > 0, use this section to 
filter the data. If not, continue 
with the program. 

Get the average of each set of 
counts using the information in 
the Bin_u arrays. 

First, set the average=0. 


! Get a count and convert it to 
! a real number. 


Use the Bin_u arrays to find out how 
many samples there are with that 
count and convert that number to a 
real number. 


Multiply the count by 

the number of times it appears in 

the data. 


! Keep a running sum of the product of 
! the count times the number of times 
! that it appears. Store this sum in 
! the place we will eventually use to 
! store the average. 
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7770 

ISOURCE 

j 



7780 

ISOURCE 


LDA 

= Ns 

7790 

ISOURCE 


STA 

Oper_l 

7800 

ISOURCE 


LDA 

=Xvar 

7810 

ISOURCE 


STA 

Resu 1 t 

7820 

ISOURCE 


JSM 

Int_to_rel 

7830 

ISOURCE 

! 



7840 

ISOURCE 


LDB 

Count 

7850 

ISOURCE 


SBL 

2 

7860 

ISOURCE 


ADB 

=Avg_u 

7870 

ISOURCE 


STB 

Oper_l 

7880 

ISOURCE 


STB 

Resul t 

7890 

ISOURCE 


LDB 

=Xvar 

7900 

ISOURCE 


STB 

Oper 2 

7910 

ISOURCE 


LDA 

= 2 

7920 

ISOURCE 


LDB 

=147155B 

7930 

ISOURCE 


JSM 

Re1_math 

7940 

ISOURCE 

! 



7950 

ISOURCE 


LDfl 

Count 

7960 

ISOURCE 


ADA 

= 1 

7970 

ISOURCE 


STA 

Count 

7980 

ISOURCE 


ADA 

=-3 

7990 

ISOURCE 


SZA 

*+2 

8000 

ISOURCE 


JMP 

Get_avg 

8010 

ISOURCE 

! 



8020 

ISOURCE 


LDA 

= 0 

8030 

ISOURCE 


STA 

Count 

8040 

ISOURCE 

Get_sdevS 

LDA 

=Zero 

8050 

ISOURCE 


LDB 

Count 

8060 

ISOURCE 


SBL 

2 

8070 

ISOURCE 


STB 

Offset 

8080 

ISOURCE 


ADB 

=De v u 

8090 

ISOURCE 


XFR 

4 

8100 

ISOURCE 


LDA 

Offset 

8110 

ISOURCE 


ADA 

= Avg_u 

8120 

ISOURCE 


STA 

I nt 

8130 

ISOURCE 


LDA 

Count 

8140 

ISOURCE 


ADA 

= 1 

8150 

ISOURCE 


LDB 

= 1024 

8160 

ISOURCE 


MPY 


8170 

ISOURCE 


ADA 

=Bi n_u 

8180 

ISOURCE 


STA 

Address 

8190 

ISOURCE 


LDA 

= 1023 

8200 

ISOURCE 


STA 

I 

8210 

ISOURCE 

! 



8220 

ISOURCE 

Dev_count : 

LDA 

= 1 

8230 

ISOURCE 


STA 

Oper_l 

8240 

ISOURCE 


LDA 

=Xvar 

8250 

ISOURCE 


STA 

Resu 1 t 

8260 

ISOURCE 


JSM 

Int_to_rel 

8270 

ISOURCE 

j 



8280 

ISOURCE 


LDA 

=Xvar 

8290 

ISOURCE 


STA 

Oper_l 

8300 

ISOURCE 


LDA 

I nt 

8310 

ISOURCE 


STA 

0per_2 

8320 

ISOURCE 


LDA 

=Yvar 

8330 

ISOURCE 


STA 

Resul t 

8340 

ISOURCE 


LDA 

= 2 

8350 

ISOURCE 


LDB 

= 1 4671 7B 

8360 

ISOURCE 


JSM 

Rel math 


! Convert the number of samples to 
! a real number. 


Divide the sum stored in Avg_u by 
the number of samples to get the 
average value of the counts. 


Get the standard deviation of each 
set of counts using the average and 
the information in the Bin_u arrays. 
First, set the standard deviation =0 


! Get a count and convert it to a 
! real number. 


! Subtract the average count from the 
! count we just got. 
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8370 

ISOURCE ! 



8380 

ISOURCE 

LDfl 

=Yvar 

8390 

ISOURCE 

STfl 

Oper_l 

8400 

ISOURCE 

STR 

0per_2 

8410 

ISOURCE 

LDfl 

=Xvar 

8420 

ISOURCE 

STR 

Resu 1 t 

8430 

ISOURCE 

LDfl 

= 2 

8440 

ISOURCE 

LDB 

= 1 47037B 

8450 

ISOURCE 

JSM 

Rel_math 

8460 

ISOURCE .! 



8470 

ISOURCE 

LDfl 

Address 

8480 

ISOURCE 

STfl 

Oper_l 

8490 

ISOURCE 

LDfl 

=Y var 

8500 

ISOURCE 

STfl 

Resul t 

8510 

ISOURCE 

JSM 

Int_to_rel 

8520 

ISOURCE ! 



8530 

ISOURCE 

LDfl 

=Xvar 

8540 

ISOURCE 

STfl 

Oper_l 

8550 

ISOURCE 

LDfl 

=Yvar 

8560 

ISOURCE 

STfl 

0per_2 

8570 

ISOURCE 

STfl 

Resu 1 t 

8580 

ISOURCE 

LDfl 

= 2 

8590 

ISOURCE 

LDB 

=147037B 

8600 

ISOURCE 

JSM 

Re1_math 

8610 

ISOURCE ! 



8620 

ISOURCE 

LDfl. 

=Yvar 

8630 

ISOURCE 

STfl 

Oper_l 

8640 

ISOURCE 

LDfl 

=Dev_u 

8650 

ISOURCE 

ADA 

Offset 

8660 

ISOURCE 

STfl 

0per_2 

8670 

ISOURCE 

STfl 

Resul t 

8680 

ISOURCE 

LDfl 

=2 

8690 

ISOURCE 

LDB 

=146721B 

8700 

ISOURCE 

JSM 

Re1_math 

8710 

ISOURCE ! 



8720 

ISOURCE 

DSZ 

Address 

8730 

ISOURCE 

DSZ 

I 

8740 

ISOURCE 

JMP 

Dev_count 

8750 

ISOURCE ! 



8760 

ISOURCE 

LDfl 

=Ns 

8770 

ISOURCE 

STfl 

Oper_l 

8780 

ISOURCE 

LDfl 

=Xvar 

8790 

ISOURCE 

STfl 

Resu 1 t 

8800 

ISOURCE 

JSM 

Int_to_rel 

8810 

ISOURCE ! 



8820 

ISOURCE 

LDB 

Offset 

8830 

ISOURCE 

flDB 

=Dev_u 

8840 

ISOURCE 

STB 

Oper_l 

8850 

ISOURCE 

STB 

Resul t 

8860 

ISOURCE 

LDB 

=Xvar 

8870 

ISOURCE 

STB 

Oper 2 

8880 

ISOURCE 

LDfl 

= 2 

8890 

ISOURCE 

LDB 

= 1 471 55B 

8900 

ISOURCE 

JSM 

Rel_math 

8910 

ISOURCE ! 



8920 

ISOURCE 

LDfl 

= Dev u 

8930 

ISOURCE 

ADA 

Offset 

8940 

ISOURCE 

STfl 

Oper_l 

8950 

ISOURCE 

STfl 

Resul t 

8960 

ISOURCE 

LDfl 

= 1 


Square the difference between the 
average and the count we got to 
get the deviation from the mean. 


Use the Bin_u arrays to find out 
how many counts there are with the 
value we have chosen, and convert 
that number to a real number. 


! Multiply the number of counts by 
I the deviation from the mean. 


Add the product to Dev_u, which 
currently contains a running sum of 
the number of counts times the 
deviation from the mean of each 
count . 


! Convert the number of samples to 
! a real number. 


Divide the sum stored in Dev_u by 
the number of samples to get the 
variance of the counts. 


Take the square root of the 
variance to get the standard 
dev i at i on . 
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8970 

ISOURCE 


LDB 

=31450B 

8980 

ISOURCE 


JSM 

Rel_math 

8990 

ISOURCE 

! 



9000 

ISOURCE 


LDA 

Count 

9010 

ISOURCE 


ADA 

= 1 

9020 

ISOURCE 


STA 

Count 

9030 

ISOURCE 


ADA 

=-3 

9040 

ISOURCE 


S2A 

*+2 

9050 

ISOURCE 


JMP 

Get_sdev 

9060 

ISOURCE 

! 



9070 

ISOURCE 


LDA 

= 0 

9080 

ISOURCE 


STA 

Count 

9090 

ISOURCE 

Dev_f i 1 tr: 

LDA 

Count 

9100 

ISOURCE 


ADA 

= 1 

9110 

ISOURCE 


LDB 

= 1024 

9120 

ISOURCE 


MPY 


9130 

ISOURCE 


ADA 

=Bi n_u 

9140 

ISOURCE 


STA 

Address 

9150 

ISOURCE 


LDA 

= 1023 

9160 

ISOURCE 


STA 

I 

9170 

ISOURCE 

! 



9180 

ISOURCE 


LDA 

Count 

9190 

ISOURCE 


SAL 

2 

9200 

ISOURCE 


ADA 

=Dev_u 

9210 

ISOURCE 


STA 

Oper_l 

9220 

ISOURCE 


LDA 

=Sdev 

9230 

ISOURCE 


STA 

0per_2 

9240 

ISOURCE 


LDA 

=Xvar 

9250 

ISOURCE 


STA 

Resul t 

9260 

ISOURCE 


LDA 

=2 

9270 

ISOURCE 


LDB 

=147037B 

9280 

ISOURCE 


JSM 

Rel _math 

9290 

ISOURCE 

1 



9300 

ISOURCE 


LDA 

=Xvar 

9310 

ISOURCE 


STA 

Qper_l 

9320 

ISOURCE 


LDA 

= Int 

9330 

ISOURCE 


STA 

Resul t 

9340 

ISOURCE 


JSM 

Re1_to_i nt 

9350 

ISOURCE 

! 



9360 

ISOURCE 


LDA 

Count 

9370 

ISOURCE 


SAL 

2 

9380 

ISOURCE 


ADA 

=Avg_u 

9390 

ISOURCE 


STA 

Oper_l 

9400 

ISOURCE 


LDA 

=Off set 

9410 

ISOURCE 


STA 

Resul t 

9420 

ISOURCE 


JSM 

Rel_to_i nt 

9430 

ISOURCE 

! 



9440 

ISOURCE 


LDB 

Offset 

9450 

ISOURCE 


TCB 


9460 

ISOURCE 

F i 1 1 r_dev: 

LDA 

I 

9470 

ISOURCE 


ADA 

B 

9480 

ISOURCE 


SAM 

*+2 

9490 

ISOURCE 


TCA 


9500 

ISOURCE 


ADA 

Int 

9510 

ISOURCE 


SAP 

*+3 

9520 

ISOURCE 


LDA 

Address 

9530 

ISOURCE 


CLR 

1 

9540 

ISOURCE 


DS2 

Address 

9550 

ISOURCE 


DS2 

I 

9560 

ISOURCE 


JMP 

F i 1 1 r dev 


Now use the standard deviation to 
filter out all the counts whose 
value is more than Sdev standard 
deviations away from the mean. 


58 



9570 

ISOURCE 

1 




9580 

ISOURCE 


I sz 

Count 

! Go through this whole process 

9590 

ISOURCE 


LDfi 

=-3 

! of checking the deviation and 

9600 

ISOURCE 


flDfi 

Count 

! filtering the data three times, 

9610 

ISOURCE 


SZfl 

*+2 

! once for each channel. 

9620 

ISOURCE 


JMP 

Dev_f i 1 tr . 


9630 

ISOURCE 

j 




9640 

ISOURCE 

j 

The 1 oop Begi 

n is performed three times, once for U, V, 

9650 

ISOURCE 

! 

and W, Each time through 

the data array is read, a count <a raw 

9660 

ISOURCE 

i 

datum) is taken from it 

and converted into a velocity. The vel- 

9670 

ISOURCE 

! 

ocity is stored in U, V, 

or M depending on whether this is the 

9680 

ISOURCE 

j 

first;, seconc 

, or third 

iteration of the loop. The first time a 

9690 

ISOURCE 

j 

particular count is encountered, the velocity correspondi ng to 

9700 

ISOURCE 

! 

it is calculated using the intermediate values found in Get_i nt , 

9710 

ISOURCE 

! 

and the velocity is stored in a table. If that count is found 

9720 

ISOURCE 

! 

again in the 

data array, 

the c orrespond i ng velocity is looked 

9730 

ISOURCE 

! 

up rather than being calculated again. 

9740 

ISOURCE 

j 




9750 

ISOURCE 

Set begin: LDfi 

= 0 


9760 

ISOURCE 


STfi 

Count 


9770 

ISOURCE 


LDfi 

Ns 


9780 

ISOURCE 


STfi 

I 


9790 

ISOURCE 


STfi 

Ngs 


9800 

ISOURCE 

! 




9810 

ISOURCE 

Begin: LDfi 

Count 


9820 

ISOURCE 


STfl 

B 


9830 

ISOURCE 


SBL 

2 


9840 

ISOURCE 


STB 

Offset 


9850 

ISOURCE 


LDB 

Ns 

! . Figure out which element of the 

9860 

ISOURCE 


MPY 


! data array we want to pick up. 

9870 

ISOURCE 


fiDfl 

I 


9880 

ISOURCE 


ADA 

=-l 


9890 

ISOURCE 


STfl 

El ementd 


9900 

ISOURCE 


LDfi 

= I nt 

! Get a raw datum from the data array. 

9910 

ISOURCE 


LDB 

=flrrayd 


9920 

ISOURCE 


JSM 

Get_e 1 ement 


9930 

ISOURCE 


LDfi 

Int 


9940 

ISOURCE 


LDB 

= 1023 

! Strip off the first six 

9950 

ISOURCE 


AND 

B 

! bits of the raw data word. 

9960 

ISOURCE 


LDB 

Count 

! See if Count = Check. 

9970 

ISOURCE 


TCB 



9980 

ISOURCE 


flDB 

Check 


9990 

ISOURCE 


SZB 

Strai ght 

! If true, use the modified data 

10000 

ISOURCE 


TCfi 


! word as an index. If not, use 

10010 

ISOURCE 


ADR 

= 1024 

! 1024 minus the data word. 

10020 

ISOURCE 

Straight: STfl 

Int 

! Store the count we have gotten. 

10030 

ISOURCE 

! 




10040 

ISOURCE 


LDfi 

Count 

! Now look in the appropriate part 

10050 

ISOURCE 


LDB 

= 1024 

! of the arrays created by the filter 

10060 

ISOURCE 


MPY 


! section. 

10070 

ISOURCE 


ADR 

Int 


10080 

ISOURCE 


fiDfl 

=-l 


10090 

ISOURCE 


ADA 

=B i n u 


10100 

ISOURCE 


LDB 

A, I 


10110 

ISOURCE 


RZB 

Goode ount 

! If the word corresponding to the 

10120 

ISOURCE 


STB 

Count 

! count is zero, the count was 

10130 

ISOURCE 


DSZ 

Ngs 

! filtered out and should be ignored. 

10140 

ISOURCE 


DSZ 

I 

! Ignore these three counts, go back 

10150 

ISOURCE 


JMP 

Beg i n 

! and do the next three. 

10160 

ISOURCE 

j 
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10170 

ISOURCE 

Goode ount : 

LDFI 

Count 

10180 

ISOURCE 


LDB 

=4896 

10190 

ISOURCE 


MPY 


10200 

ISOURCE 


LDB 

Int 

10210 

ISOURCE 


FlDB 

= -l 

10220 

ISOURCE 


SBL 

2 

10230 

ISOURCE 


RDF) 

B 

10240 

ISOURCE 


RDfl 

= 1 

10250 

ISOURCE 


ADA 

=flrray 1 

10260 

ISOURCE 


STR 

Address 

10270 

ISOURCE 


LDB 

= I nt 2 

10280 

ISOURCE 


XFR 

1 

10290 

ISOURCE 


LDFI 

Int2 

10300 

ISOURCE 


SZfl 

Cal cul ate 

10310 

ISOURCE 


JMP 

Over 

10320 

ISOURCE 

Cal cul ate: 

LDFI 

= I nt 

10330 

ISOURCE 


STfl 

Oper_l 

10340 

ISOURCE 


LDfi 

=Y var 

10350 

ISOURCE 


STR 

Resul t 

10360 

ISOURCE 


JSM 

Int_t o_rel 

10370 

ISOURCE 


STB 

Oper 2 

10380 

ISOURCE 


LDfl 

=R1 

10390 

ISOURCE 


RDfl 

Offset 

10400 

ISOURCE 


STR 

Oper_l 

10410 

ISOURCE 


LDfl 

=Xvar 

10420 

ISOURCE 


STR 

Resul t 

10430 

ISOURCE 


LDfl 

=2 

10440 

ISOURCE 


LDB 

=147155B 

10450 

ISOURCE 


JSM 

Rel_math 

10460 

ISOURCE 

! 



10470 

ISOURCE 


LDfl 

=Numi xl 

10480 

ISOURCE 


ADA 

Offset 

18490 

ISOURCE 


STfl 

Oper_l 

10500 

ISOURCE 


LDfl 

=Xvar 

10510 

ISOURCE 


STfl 

0per_2 

10520 

ISOURCE 


LDfl 

=Yvar 

10530 

ISOURCE 


STfl 

Resu 1 t 

10540 

ISOURCE 


LDfl 

=2 

10550 

ISOURCE 


LDB 

=146717B 

10560 

ISOURCE 


JSM 

Rel_math 

10570 

ISOURCE 

! 



10580 

ISOURCE 


LDB 

Count 

10590 

ISOURCE 


flDB 

= -l 

10600 

ISOURCE 


SBM 

*+2 

10610 

ISOURCE 


JMP 

Sames i gn 

10620 

ISOURCE 


LDfl 

=2ero 

10630 

ISOURCE 


STfl 

Oper_l 

10640 

ISOURCE 


LDfl 

Resul t 

10650 

ISOURCE 


STfl 

0per_2 

10660 

ISOURCE 


LDfl 

=Xvar 

10670 

ISOURCE 


STfl 

Resul t 

10680 

ISOURCE 


LDfl 

=2 

10690 

ISOURCE 


LDB 

= 1 4671 7B 

10700 

ISOURCE 


JSM 

Rel_math 

10710 

ISOURCE 

Sames i gn: 

NOP 


10720 

ISOURCE 

! 



10730 

ISOURCE 


LDfl 

Resul t 

10740 

ISOURCE 


STfl 

Oper 1 

10750 

ISOURCE 


LDfl 

= Df 1 

10760 

ISOURCE 


ADA 

Of fs=* 


Now, use Count to find out 
which lookup table array 
we want to use, and use the 
count we got from the data array 
to find exactly where in the 
table we want to go. 


! If that table entry is zero, 
! calculate a velocity for it. 


! Convert the count into 
! a real number. 


Divide the range we found 
earlier by the count to get a 
frequency. 


! Find < Nub-Num i xN ) -Frequenc yN . 


If we are calculating U, reverse 
the sign of < (Nub-Numi xN)-FrequencyN> 
so as to reverse the sign of U. 

Leave V and Wv alone. 


Find Ve 1 oc i t y= ( < Nub-Num i xN > 
-FrequencyH>*<Mi 1 1 *DfN> 
and store in a place in the 
lookup table c orrespond i ng to 
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10770 

ISOURCE 


STfi 

0per_2 ! 

the data count. 

10780 

ISOURCE 


LDfl 

Address 


10790 

ISOURCE 


ADR 

=-l 


10800 

ISOURCE 


STR 

Resul t 


10810 

ISOURCE 


LDR 

=2 


10820 

ISOURCE 


LDB 

=147037B 


10830 

ISOURCE 


JSM 

Rel_math 


1084O 

ISOURCE 

! 




10350 

ISOURCE 

Over: LDR 

Rddress ! 

Transfer the velocity from 

10860 

ISOURCE 


RDR 

=-l ! 

the lookup table to U, V, or W, 

10870 

ISOURCE 


LDB 

=U ! 

as appropriate. 

10880 

ISOURCE 


RDB 

Offset 


10890 

ISOURCE 


XFR 

4 


10900 

ISOURCE 

! 




10910 

ISOURCE 


ISZ 

Count 


10920 

ISOURCE 


LDfl 

= t3 ! 

Have U, V, and W al 1 

10930 

ISOURCE 


RDR 

Count ! 

been calculated? 

10940 

ISOURCE 


SZR 

*+2 ! 

If not, go back again. 

10950 

ISOURCE 


JMP 

Begi n 


10960 

ISOURCE 


STR 

Count ! 

If so, set Count = 0. 

10970 

ISOURCE 

i 




10980 

ISOURCE 

i 

Now we convert the W we have obtained (which is measured at 

10990 

ISOURCE 

I 

ang 1 e Theta 

to the V-axis) 

to the W we want (which should 

11000 

ISOURCE 

j 

be measured 

at an angle of 

90 degrees to the V-axis). 

11010 

ISOURCE 

j 

Thus find W= 

<Wv-V*Cos<Theta))/'Sin<Theta>. 

11020 

ISOURCE 

j 




11030 

ISOURCE 


LDR 

=V 


11040 

ISOURCE 


STR 

Oper_l 


11050 

ISOURCE 


LDR 

=Cos 


11060 

ISOURCE 


STR 

0per_2 


11070 

ISOURCE 


LDfl 

= Xvar 


11080 

ISOURCE 


STR 

Resul t 


11090 

ISOURCE 


LDR 

= 2 


1 1 100 

ISOURCE 


LDB 

= 1 47037B 


11110 

ISOURCE 


JSM 

Rel_math 


11120 

ISOURCE 


LDfl 

Resul t 


11130 

ISOURCE 


STR 

Oper 2 


1 1 140 

ISOURCE 


LDfl 

=W 


11150 

ISOURCE 


STR 

Oper_l 


11160 

ISOURCE 


LDR 

=Yvar 


11170 

ISOURCE 


STR 

Resul t 


1 1 180 

ISOURCE 


LDfl 

=2 


1 1 190 

ISOURCE 


LDB 

=146717B 


11200 

ISOURCE 


JSM 

Rel_math 


11210 

ISOURCE 


LDfl 

Resul t 


11220 

ISOURCE 


STR 

Oper_l 


11230 

ISOURCE 


LDR 

=Si n 


11240 

ISOURCE 


STR 

Oper 2 


11250 

ISOURCE 


LDR 

=W 


11260 

ISOURCE 


STR 

Resul t 


11270 

ISOURCE 


LDR 

= 2 


11280 

ISOURCE 


LDB 

=147155B 


11290 

ISOURCE 


JSM 

Rel_math 


11300 

ISOURCE 

j 




11310 

ISOURCE 

! 

Now take running sums of U, 

V, kl, and several products 

11320 

ISOURCE 

j 

of these velocities. The sums are taken using the utility 

11330 

ISOURCE 

j 

"Add". The sums are calculated in an unusual sequence in 

11340 

ISOURCE 

i 

order to reduce the number 

of program steps needed to cal- 

11350 

ISOURCE 

! 

culate them. 



11360 

ISOURCE 

j 





an 
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11370 ISOURCE 
11380 ISOURCE 
11390 ISOURCE 
11400 ISOURCE ! 
11410 ISOURCE 
11420 ISOURCE 
11430 ISOURCE 
11440 ISOURCE 
11450 ISOURCE 
11460 ISOURCE 
11470 ISOURCE 
11480 ISOURCE ! 
11490 ISOURCE 
11500 ISOURCE 
11510 ISOURCE 
11520 ISOURCE ! 
11530 ISOURCE 
11540 ISOURCE 
11550 ISOURCE 
11560 ISOURCE 
11570 ISOURCE 
11580 ISOURCE 
11590 ISOURCE 
11600 ISOURCE ! 
11610 ISOURCE 
11620 ISOURCE 
11630 ISOURCE 
11640 ISOURCE ! 
11650 ISOURCE 
11660 ISOURCE 
11670 ISOURCE 
11680 ISOURCE 
11690 ISOURCE 
11780 ISOURCE 
11710 ISOURCE 
11720 ISOURCE 
11730 ISOURCE 
11740 ISOURCE ! 
11750 ISOURCE 
11760 ISOURCE 
11770 ISOURCE 
11780 ISOURCE ! 
11790 ISOURCE 
11800 ISOURCE 
11810 ISOURCE 
11820 ISOURCE ! 
11830 ISOURCE 
11840 ISOURCE 
11850 ISOURCE 
11860 ISOURCE 
11870 ISOURCE 
11880 ISOURCE 
11890 ISOURCE 
11900 ISOURCE ! 
11910 ISOURCE 
11920 ISOURCE 
11930 ISOURCE 
11940 ISOURCE ! 
11950 ISOURCE 
11960 ISOURCE 


LDfl 

=Su 

LDB 

=U 

JSM 

Add 

LDfl 

Qpen_2 

STfl 

Open 1 

LDfl 

=Uu 

STfl 

Resul t 

LDfl 

=2 

LDB 

=147037B 

JSM 

Re1_math 

LDfl 

=Suu 

LDB 

Resul t 

JSM 

Add 

LDfl 

= V 

STfl 

Oper_l 

LDfl 

= Xuan 

STfl 

Resul t 

LDfl 

= 2 

LDB 

=147037B 

JSM 

Rel_math 

LDfl 

=Suuu 

LDB 

Resul t 

JSM 

Add 

LDfl 

=Uu 

STfl 

Open 1 

LDfl 

=W 

STfl 

0pen_2 

LDfl 

=Xuan 

STfl 

Resu 1 t 

LDfl 

=2 

LDB 

= 1 47037B 

JSM 

Re1_math 

LDfl 

= Suuw 

LDB 

Resul t 

JSM 

Add 

LDfl 

=Sv 

LDB 

=V 

JSM 

Add 

LDfl 

=U 

STfl 

Open 1 

LDfl 

=Uu 

STfl 

Resul t 

LDfl 

=2 

LDB 

=1470373 

JSM 

Rel_math 

LDfl 

=Suv 

LDB 

Resul t 

JSM 

Add 

LDfl 

=V 

STfl 

Open_l 


! Find Su=Su+U 


! Find U*U 


! Find Suu=Suu+<U*U> 


! Find U*U*V 


! Find Suuu=Suuu+<U*U*V:> 


! Find U*U*W 


! Find Suuu=Suuw+ < U*U*W > 


! Find Su = Sv+V 


! Find U*V 


! Find Suu=Suy+(U*V> 


! Find U*V*V 
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11970 

11980 

11990 

12000 

12010 

12020 

12030 

12040 

12050 

12060 

12070 

12080 

12090 

12100 

12110 

12120 

12130 

12140 

12150 

12160 

12170 

12180 

12190 

12200 

12210 

12220 

12230 

12240 

12250 

12260 

12270 

12280 

12290 

12300 

12310 

12320 

12330 

12340 

12350 

12360 

12370 

12380 

12390 

12400 

12410 

12420 

12430 

12440 

12450 

12460 

12470 

12480 

12490 

12500 

12510 

12520 

12530 

12540 

12550 

12560 


ISOURCE 

LDfl 

=Xuar 

I SOURCE 

STfl 

Resu 1 t 

ISOURCE 

LDfl 

=2 

ISOURCE 

LDB 

= 1 47037B 

ISOURCE 
ISOURCE ! 

JSM 

Rel math 

ISOURCE 

LDfl 

=Suuu 

ISOURCE 

LDB 

Resul t 

ISOURCE 
ISOURCE ! 

JSM 

Add 

ISOURCE 

LDfl 

=Uu 

ISOURCE 

STfl 

Open 1 

ISOURCE 

LDfl 

=W 

ISOURCE 

STfl 

0per_2 

ISOURCE 

LDfl 

=Xuar 

ISOURCE 

STfl 

Resul t 

ISOURCE 

LDfl 

=2 

ISOURCE 

LDB 

=147037B 

ISOURCE 
ISOURCE ! 

JSM 

Rel_math 

ISOURCE 

LDfl 

=Suuw 

ISOURCE 

LDB 

Resul t 

ISOURCE 

JSM 

Add 

ISOURCE ! 



ISOURCE 

LDfl 

=V 

ISOURCE 

STfl 

Oper_l 

ISOURCE 

STfl 

0per_2 

ISOURCE 

LDfl 

=Xuar 

ISOURCE 

STfl 

Resul t 

ISOURCE 

LDfl 

=2 

ISOURCE 

LDB 

= 1 47037B 

ISOURCE 
ISOURCE ! 

JSM 

Rel_math 

ISOURCE 

LDfl 

=Suu 

ISOURCE 

LDB 

Resul t 

ISOURCE 
ISOURCE ! 

JSM 

Add 

ISOURCE 

LDfl 

=Sw 

ISOURCE 

LDB 

=W 

ISOURCE 
ISOURCE ! 

JSM 

Add 

ISOURCE 

LDfl 

=U 

ISOURCE 

STfl 

Oper_l 

ISOURCE 

LDfl 

=Xuar 

ISOURCE 

STfl 

Resul t 

ISOURCE 

LDfl 

=2 

ISOURCE 

LDB 

=147037B 

ISOURCE 
ISOURCE ! 

JSM 

Rel_math 

ISOURCE 

LDfl 

=Suw 

ISOURCE 

LDB 

Resu 1 t 

ISOURCE 
ISOURCE ! 

JSM 

Add 

ISOURCE 

LDfl 

=W 

ISOURCE 

STfl 

Open 1 

ISOURCE 

LDfl 

=V 

ISOURCE 

STfl 

0per_2 

ISOURCE 

LDfl 

=Xuar 

ISOURCE 

STfl 

Resul t 

ISOURCE 

LDfl 

=2 


! Find Suvv=Suyy+<U*V*V) 


! Find U*V*W 


! Find Suvw=Suuw+<U*V*W) 


! Find V*V 


! Find Sw=Sw+<V*V) 


! Find Sw=Sw+W 


! Find U*W 


! Find Suw=Suw+<U*W) 


! Find V*W 
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12570 

ISOURCE 



LDB 

= 1 47037B 



12580 

1 SOURCE 



JSM 

Re 1 _math 



12590 

ISOURCE 

j 






12600 

ISOURCE 



LDfl 

=Svw 

! Find 

Svw=Svw+(V*W) 

12610 

ISOURCE 



LDB 

Resul t 



12620 

ISOURCE 



JSM 

Add 



12630 

ISOURCE 

j 






12640 

ISOURCE 



LDfl 

=W 

! Find 

w*w 

12650 

ISOURCE 



STfl 

Oper_l 



12660 

ISOURCE 



STfi 

0per_2 



12670 

ISOURCE 



LDfl 

=Xvar 



12680 

ISOURCE 



STfl 

Resul t 



12690 

ISOURCE 



LDfl 

=2 



12700 

ISOURCE 



LDB 

= 1 47037B 



12710 

ISOURCE 



JSM 

Rel_math 



12720 

ISOURCE 

! 






12730 

ISOURCE 



LDfl 

=Sww 

! Find 

Sww=Sww+<W*W) 

12740 

ISOURCE 



LDB 

Result 



12750 

ISOURCE 



JSM 

Add 



12760 

ISOURCE 

j 






12770 

ISOURCE 



LDfl 

=U 

! Find 

U*W*W 

12780 

ISOURCE 



STfl 

Oper 1 



12790 

ISOURCE 



LDfl 

=Uu 



12800 

ISOURCE 



STfl 

Resul t 



12810 

ISOURCE 



LDfl 

=2 



12820 

ISOURCE 



LDB 

=147037B 



12830 

ISOURCE 



JSM 

Re1_math 



12840 

ISOURCE 

! 






12850 

ISOURCE 



LDfl 

=Suww 

! Find 

Suww=Suww+<U*W*U> 

12860 

ISOURCE 



LDB 

Resul t 



12870 

ISOURCE 



JSM 

Add 



12880 

ISOURCE 

! 






12890 

ISOURCE 



DS2 

I 

! Continue to calculate running 

12900 

ISOURCE 



JMP 

Begi n 

! sums 

until out of samples. 

12910 

ISOURCE 

j 






12920 

ISOURCE 

j 

Now 

place the finished 

summs i n 

the COMMON region so that 

12930 

ISOURCE 

j 

the 

BASIC program has 

access to 

them, and then return to 

12940 

ISOURCE 

j 

the 

BASIC program. 



12950 

ISOURCE 

Replace; LDfl 

=Rangel 



12960 

ISOURCE 



LDB 

=Ran l_par 



12970 

ISOURCE 



JSM 

Put_val ue 



12980 

ISOURCE 



LDfl 

=Range2 



12990 

ISOURCE 



LDB 

=Ran2_par 



13000 

ISOURCE 



JSM 

Put_val ue 



13010 

ISOURCE 



LDfl 

=Range3 



13020 

ISOURCE 



LDB 

=Ran3_par 



13030 

ISOURCE 



JSM 

Put value 



13040 

ISOURCE 



LDfl 

=Su 



13050 

ISOURCE 



LDB 

=Su_par 



13060 

ISOURCE 



JSM 

Put value 



13070 

ISOURCE 



LDfl 

=S v 



13080 

ISOURCE 



LDB 

=Sv_par 



13090 

ISOURCE 



JSM 

Put_va1 ue 



13100 

ISOURCE 



LDfl 

=Sw 



13110 

ISOURCE 



LDB 

=Sw_par 



13120 

ISOURCE 



JSM 

Put_val ue 



13130 

ISOURCE 



LDfl 

=Suu 



13140 

ISOURCE 



LDB 

=Suu_par 



13150 

ISOURCE 



JSM 

Put_val ue 



13160 

ISOURCE 



LDfl 

=Svv 
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13170 

13180 

13190 

13200 

13210 

13220 

13230 

13240 

13250 

13260 

13270 

13280 

13290 

13300 

13310 

13320 

13330 

13340 

13350 

13360 

13370 

13380 

13390 

13400 

13410 

13420 

13430 

13440 

13450 

13460 

13470 

13480 

13490 

13500 

13510 

13520 

13530 

13540 

13550 

13560 

13570 

13580 

13590 

13600 

13610 

13620 

13630 

13640 

13650 

13660 

13670 

13680 

13690 

13700 

13710 

13720 

13730 

13740 

13750 

13760 


ISOURCE 


LDB 

=Suu_par 

ISOURCE 


JSM 

Put_ual ue 

ISOURCE 


LDfl 

=Sww 

ISOURCE 


LDB 

=Swu_par 

ISOURCE 


JSM 

Put_ual ue 

ISOURCE 


LDfl 

= Suu 

ISOURCE 


LDB 

=Suv_par 

ISOURCE 


JSM 

Put_ual ue 

ISOURCE 


LDfl 

=Suw 

ISOURCE 


LDB 

=Suw_par 

ISOURCE 


JSM 

Put_ual ue 

ISOURCE 


LDfl 

=Suw 

ISOURCE 


LDB 

=Sww_par 

ISOURCE 


JSM 

Put_ual ue 

ISOURCE 


LDfl 

=Suuu 

ISOURCE 


LDB 

=Suuw_par 

ISOURCE 


JSM 

Put_ual ue 

ISOURCE 


LDfl 

= Suuu 

ISOURCE 


LDB 

=Suwu_par 

ISOURCE 


JSM 

Put_wal ue 

ISOURCE 


LDfl 

= Suuw 

ISOURCE 


LDB 

=Suuw_par 

ISOURCE 


JSM 

Put_ual ue 

ISOURCE 


LDfl 

=Suww 

ISOURCE 


LDB 

=Suww_par 

ISOURCE 


JSM 

Put_wal ue 

ISOURCE 


LDfl 

=Suuu 

ISOURCE 


LDB 

=Suvw_par 

ISOURCE 


JSM 

Put_ual ue 

ISOURCE 


LDfl 

=Ngs 

ISOURCE 


LDB 

=Ngs_par 

ISOURCE 

ISOURCE 

! 

JSM 

Put_ual ue 

ISOURCE 


LDfl 

= B i n u 

ISOURCE 


LDB 

= 192 

ISOURCE 


STB 

I 

ISOURCE 


LDB 

= 16 

ISOURCE 


STB 

Count 

ISOURCE 


LDB 

Cbin par 

ISOURCE 

Transf er: 

XFR 

16 

ISOURCE 


ADA 

Count 

ISOURCE 


ADB 

Count 

ISOURCE 


DSZ 

I 

ISOURCE 


JMP 

Transf er 

ISOURCE 

ISOURCE 

i 

RET 

1 

ISOURCE 

ISOURCE 

! 

LIT 

200 

ISOURCE 

ISOURCE 

! The u> i 1 
! 

i ty 

'Add" is us 

ISOURCE 

Add: 

ISZ 

Ut 1 count 

ISOURCE 


STfl 

Oper_l 

ISOURCE 


STB 

0per_2 

ISOURCE 


STfl 

Resul t 

ISOURCE 


LDfl 

= 2 

ISOURCE 


LDB 

= 1 4672 1 B 

ISOURCE 


JSM 

Rel_math 

ISOURCE 


DSZ 

Ut 1 count 

ISOURCE 


RET 

1 

ISOURCE 


JSM 

Ut 1 end 


Transfer the contents of Bin_n 
to the common area without using 
the slow HP-supplied external 
subrout i nes . 


to add up the running sums. 
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13770 ISOURCE ! 

13780 ISOURCE END Find vel 
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TABLE 1 

MACRODYNE RANGE SETTINGS 

MANTISA 1 500 1024 


RANGE 


0 

32 GHz 

1 

16 " 

2 

8 " 

3 

4 " 

4 

2 " 

5 

1 ” 

6 

500 MHz 

7 

250 " 

8 

125 " 

9 

62.5 ” 

10 

31.25 ” 

11 

15.62 ” 

12 

oo 
>— >■ 

13 

3.90 ” 

14 

1.95 " 


64 MHz 31.28 MHz 


32 M 

15.64 

16 M 

7.20 

8 ” 

3.91 

4 " 

1.95 

2 ” 

977 1 

1 M 

488 

500 kHz 

244 

250 M 

122 

125 M 

61 

62.5 ” 

30 

31.25 " 

15 

15.62 ” 

7.6 

7.81 " 

3.8 

3.90 " 

1.9 
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4-WATT 

ARGON-ION 

LASER 


2 COLOR FOUR 




BEAM MATRIX f \ 

TO TRANSMITTING \ 7 

TABLE 


POLARIZATION 
ROTATOR 


( 2 ) 


m 


ROTATING 

PRISM 

BEAM 

DISPLACER 

BEAM 

STEERING 

MODULES 



BRAGG CELL 
(BLUE) 

BRAGG CELL 
(GREEN) 

BEAM SPLITTER 
(BLUE) 

POLARIZATION 

ROTATOR 

(BLUE) 

BEAM DISPLACER 
(BLUE) 

BEAM SPLITTER 
(GREEN) 

POLARIZATION 
ROTATOR 
(GREEN) 


(17) 


(16) 


(15) 


(9) 


( 8 ) 


(14) 


(13) 


( 12 ) 


(ID 


( 10 ) 


GREEN 514,5 nm4 i BLUE 488 nm 


MIRROR 



ATTENUATOR 

(3) 



DISPERSION 

PRISM 

(4) 


MIRROR 
(7) 


MIRROR 

( 6 ) 


COLOR SEPARATOR BOX 


1 


Schematic of the optics table layout 





Fig. 2 Schematic of the transmitting optics. 

















Fig. 4 


Schematic of the receiving optics. 












TUBE 

Fi 9* 5 Evaluation of the effective probe length. 


HI-PASS SINE UNIVERSAL 

PHOTOMULTIPLIER AMPLIFIERS FILTERS WAVE FREQUENCY SINGLE BURST DIGITAL 

TUBES (GAIN' 37 dB) (10 MHz) GEN COUNTERS COUNTERS VOLTMETERS 



Fig. 6 


Schematic of the signal processing electronics. 
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FROM 

LDV 


Fig. 


DIGITAL 



Computer interface simplified block diagram. 


TO 

HP 9845 
CPU 







Fig. 8 


Block diagram for data acquisition program. 


76 











Fig. 9 


Block diagram for data reduction program. 
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SYNC PULSES FROM COUNTER 


DIGITAL MULTIPLEXER INPUTS 16 BIT WORDS 


J3 J4 J5 J6 J7 J8 



DIGITAL DATA 
FROM COUNTER 


J9 


OUTPUT TO 
PDP-1 1/DR-1 IB 


A/O DIGITAL OUTPUTS 12 BITS 2 S COMPLIMENT 
J15 J16 J17 J18 



CHANNEL 1 CHANNEL 2 CHANNEL 3 CHANNEL 4 


J10 


hp&e::' 


J11 J12 J13 J14 

© © © © 

ANALOG INPUTS 


OUTPUT TO 
HP9845 98032A 


O ® 

POWER POWER 
FUSE INPUT 


EVENT SYNC MODE 

i~ ii ■'» -'itxar j] 

AB AC BC ABC RANDOM 


COINCIDENCE TIME 



LDV EVENT SYNCHRONIZATION 


EVENT RATE A 



COUNTER CLOCK FREQUENCY 


IBSE 


TIME INTERVAL 
COUNTER 

_ overflow 

3 o 


LDV EVENT INPUTS 

ABC 


1 

□ 

□ 

EVENT RATE B 

c 

: 

□: 

EVENT RATE C 

□ 

□ 

□ 





TIME 

COUNTER 
RESET 
EXT INT 

i 


LDV INHIBIT OUTPUTS external reset 


© © 


GATE TIME 


INTERFACE CONTROL 
EVENT MODE 
LDV LDV A/D 


[ feair: ! | 


WORDS MULTIPLEXED 




1E=<31 


COINCIDENCE RATE 


WORD STATUS 

12 3 4 5 6 7 8 

oooooooo 


0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

OOOOOOOO OOOOOOOO 

EVENT ARM 

o 

HANDSHAKE STATUS 

HP BO* 1 1 
CILO PFLG PCTL*~- 

O o 0]g=tC~~1 


RDV GO POST BUSY 

O O O O 


MANUAL EVENT RATE 
EVENT WARNING 


o 

POWER 


NASA-ARC 

30-LDV COMPUTER INTERFACE MOD B 


I^Xf] INDICATES SELECTED OPERATING MODE 


Fig. 10 NASA LDV- A/D computer interface connections 

and settings. 
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HONEYCOMB 




f 


Fig. 11 Experimental rig. 
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LULU 



(a) Cross-wire measurements 



(b) LDV measurements 


Fig. 12 


Secondary velocity plots. 


\ 
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Fig. 13 


Comparison of X-wire and LDV measurements in a 
vortex/mixing layer interaction. 
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